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FAILURE INITIATION AND PROGRESSION IN INTERNALLY PRESSURIZED 
NON-CIRCULAR COMPOSITE CYLINDERS 

ABSTRACT 

In this study, a progressive failure analysis is used to investigate leakage in internally pressur- 
ized non-circular composite cylinders. This type of approach accounts for the localized loss of 
stiffness when material failure occurs at some location in a structure by degrading the local mate- 
rial elastic properties by a certain factor. The manner in which this degradation of material proper- 
ties takes place depends on the failure modes, which are determined by the application of a failure 
criterion. The finite-element code STAGS, which has the capability to perform progressive failure 
analysis using different degradation schemes and failure criteria, is utilized to analyze laboratory 
scale, graphite-epoxy, elliptical cylinders with quasi-isotropic, circumferentially-stiff, and axially- 
stiff material orthotropies. The results are divided into two parts. The first part shows that leakage, 
which is assumed to develop if there is material failure in every layer at some axial and circumfer- 
ential location within the cylinder, does not occur without failure of fibers. Moreover before fibers 
begin to fail, only matrix tensile failures, or matrix cracking, takes place, and at least one layer in 
all three cylinders studied remain uncracked, preventing the formation of a leakage path. That 
determination is corroborated by the use of different degradation schemes and various failure cri- 
teria. Among the degradation schemes investigated are the degradation of different engineering 
properties, the use of various degradation factors, the recursive or non-recursive degradation of 
the engineering properties, and the degradation of material properties using different computa- 
tional approaches. The failure criteria used in the analysis include the noninteractive maximum 
stress criterion and the interactive Hashin and Tsai-Wu criteria. The second part of the results 
shows that leakage occurs due to a combination of matrix tensile and compressive, fiber tensile 
and compressive, and inplane shear failure modes in all three cylinders. Leakage develops after a 
relatively low amount of fiber damage, at about the same pressure for three material orthotropies, 
and at approximately the same location. 
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Chapter 1 Introduction 


Because of aerodynamic advantages, weight, and geometric constraints, future designs of fuel 
tanks for reusable launch vehicles may call for noncircular cylinders fabricated of fiber-reinforced 
laminated composite materials. A critical issue to consider ih the design of composite tanks is fail- 
ure due to leakage of fuel through the tank wall. Typical fuels used in these vehicles are gaseous 
substances, such as hydrogen and oxygen, that are liquefied at cryogenic temperatures. The inter- 
nal pressure resulting from the use of these gases can produce enough material damage in the wall 
so that a continuous path through the layers of the composite laminate can occur. This situation is 
extremely hazardous, as it involves leakage of volatile fluids. To prevent this hazard it is essential 
that the conditions under which leakage is likely to occur be known. 

To study the potential for leakage in noncircular cylinders, a progressive failure analysis will 
be used here. Such an approach is based on the manner by which composite materials are known 
to fail. Generally, at a high enough load level, small local cracks in the material initiate the failure 
process. Due to the material softening because of the cracking, stresses redistribute and more 
cracking occurs. As the load level increases, the failed, or damaged, region becomes larger and 
stresses are redistributed further. At some level of load, fibers begin to fail at a particular location. 
Stresses are again redistributed, more failures occur, the load is increased, more fibers fail, more 
matrix cracking occurs, and the region of damage continuous to grow. At some point the material 
can no longer sustain the load and the laminate breaks into two pieces, or in the case of a pressur- 
ized cylinder, the pressure rapidly escapes and a large hole develops. The work here focuses on an 
earlier step of damage, namely, one when the cylinder is intact but leakage is beginning. This 
could occur due to excessive matrix cracking with no fiber failure, or it may involve a small 
amount of fiber failure. The primary objective of the study here is to determine at what pressure 
levels leakage occurs and what the modes of failure are when this happens. Here, leakage will be 
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assumed to occur if there is material failure in every layer at some axial and circumferential loca- 
tion within the cylinder. 

To quantitatively assess the pressure levels and the modes of failure associated with leakage in 
noncircular laminated composite cylinders, it is necessary to adopt a specific type of noncircular 
geometry, rui tiiis Si'uuy, a cylinder with an elliptical cross section will uc used. The reason for 
this choice is that past works have characterized some aspects of the response of elliptical com- 
posite cylinders and, if needed, closed-form expressions for the cross-section geometry are avail- 
able. Among the past works that described the behavior of elliptical cylinders is the study carried 
out by McMurray and Hyer in [1, 2]. In that study, using a semi-analytical approach, the effects of 
noncircular geometry, geometric nonlinearities, and different material orthotropies on the 
response and failure of internally pressurized elliptical cylinders with clamped boundary condi- 
tions were investigated. To look at the effects of noncircular geometry, a portion of the study com- 
pared the geometrically linear response of a quasi-isotropic [±45/0 /90] f elliptical cylinder, 

where 0 deg. is the axial direction, to that of a circular cylinder with the same lay-up. In the mod- 
eling of the elliptical cylinder, McMurray and Hyer used Donnell shell theory. This portion of the 
study concluded that in elliptical cross-sections the displacement response varies 'with axial and 
circumferential position, producing spatially varying displacement gradients. Through the consti- 
tutive relation, that variation is reflected in the stresses, resulting in stress concentrations at spe- 
cific axial and circumferential locations. The specific locations of those stress concentrations 
found by McMurray and Hyer can be summarized as follows: 

i. The maximum force resultants in the axial direction occur at the clamped boundaries at the ends 
of the major diameter. 

ii. The maximum force resultants in the circumferential direction occur partway between the 
clamped boundaries and the midspan of the cylinder at the ends of the minor diameter. 

iii. The maximum shear force resultants occur near the clamped boundaries part way between the 
ends of the minor diameter and the ends of the major diameter. It is important to note that due to 
the presence of the stiffnesses Z) 16 and Z) 2 g, those maximum values do not occur in every quadrant 
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of the cylinder. In fact, they occur in opposite quadrants at one clamped boundary and the reverse 
of that at the other clamped boundary. 

iii. The maximum magnitude of the bending moment resultants in the axial and circumferential 
direction and shear moment resultants occur at the clamped boundaries at the ends of the minor 
diameter. 

In investigating the influence of geometric nonlinearities on the response of elliptical cylin- 
ders, McMurray and Hyer found that if geometric nonlinearities are included in the analysis, the 
magnitude and location of the maximum displacement response are slightly different with respect 
to the geometrically linear analysis. Those nonlinear effects are found in the strain and stress 
responses as well. To determine the effects of material orthotropy, the study compared the quasi- 
isotropic cylinder with two other cylinders which had circumferentially-stiflf, [±A5l90 2 ^\i2^ 311(1 
axially-stiflf, [±45/02/90 j/ 2], lay-ups. Despite the fact these cylinders were of different wall thick- 
ness, these lay-ups were chosen because for each cylinder there was at least one layer with fibers 
in the axial direction and one layer with fibers in the circumferential direction. To do otherwise 
would be unrealistic. Using a geometrically nonlinear analysis, the study found that material 
orthotropy diminishes or amplifies the effect of elliptical geometry in some of the responses. For 
example, an axially-stiff material orthotropy amplifies the circumferential strain response of an 
elliptical cross-section, while a circumferentially-stiflf material orthotropy tends to diminish it. On 
the other hand, a circumferentially-stiflf material orthotropy accentuates the axial strain response 
of an elliptical cross section, while the axially-stiflf material orthotropy mitigates it. Making use of 
the above findings, McMurray and Hyer included failure predictions. By using the maximum 
stress and Hashin failure criteria and geometrically linear and nonlinear analyses, the failure pres- 
sure, location, and mode were predicted for each of the three cylinder lay-ups. The study found 
that both criteria predicted that the initial mode of failure was matrix tension, or matrix cracking, 
and that it occurred at the clamped boundaries in all the cylinders. The two criteria predicted that 
for the quasi-isotropic and axially-stiflf cylinders the pressure for initial failure was higher than 
that for the circumferentially-stiflf cylinder. Regarding the circumferential location and layer 
where the matrix failure occurred: Both criteria predicted that, in the first two cylinders, matrix 
cracking occurred partway between the ends of the minor and major diameters at the inner most 
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+45° layer. This cracking did not occurred in all the quadrants of the cylinder due to the presence 
of the stiffnesses D\^ and Z?26- The cracking occurred at only two opposite quadrants at each 
clamped boundary. According to both criteria, first failure in the circumferentially-stiff cylinder 

appeared at the ends of the major diameter at the outer most 90° layer. When geometrical nonlin- 
earities were included in the analysis, in general, the failure pressures were higher and circumfer- 
ential locations were somewhat different than when only geometrical linearities were used. The 
study also determined the first fiber failure pressure, location, and mode using the same two fail- 
ure criteria and analyses. In this analysis of first fiber failure, the effects of any previous matrix 
failures were ignored. For the geometrically linear case, both criteria found that the first fiber fail- 
ure was due to fiber compression in the three cylinders. Although different for each degree of 
orthotropy, the two criteria predicted similar failure pressures and locations. For the geometrically 
nonlinear case, it was found that the Hashin criterion predicted that the mode of failure was fiber 
tension for the circumferentially-stiff cylinder; however, the maximum stress criterion predicted 
that the first fiber failure was due to fiber compression. 

Though useful, since the findings in [1] give an accurate idea of where and when initial matrix 
cracking may occur in elliptical composite cylinders, the effort of predicting first fiber failure was 
somewhat limited. Because no consideration was given to the loss of local stiffness near the 
region where the first matrix crack occurred and to the consequent redistribution of stresses, 
which may cause additional matrix cracks, the first fiber failure prediction may not have been 
accurate. In other words, the study did not use a progressive failure analysis to determine the 
effects of loss of local stiffness and development of additional cracks on the first fiber failure 
pressure, location, and mode of failure. In a later work, McMurray and Hyer [3] accounted for the 
loss of stiffness due to previous matrix cracks. In that study, they used a progressive failure analy- 
sis to characterize the matrix cracking between the first matrix failure and the first fiber failure 
events, and to determine if a leakage path developed at some pressure between those two events. 
The study examined the matrix cracking scenario, predicted the first fiber failure, and determined 
if a leakage path was formed for the same internally pressurized elliptical cylinders investigated 
in [1]. In this case, however, instead of using a semi -analytical approach to approximate the geo- 
metrically nonlinear response of the cylinders, McMurray and Hyer used the finite-element code 
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STAGS (STructural Analysis of General Shells) in combination with a stand-alone pre/post pro- 
cessing program to perform the progressive failure analysis. 

A typical progressive failure methodology, and the one used in [3], uses five basic steps, as 
explained by Sleight, Knight, and Wang [4], First, a nonlinear static analysis is performed to sat- 
isfy the equilibrium of the structure which is under a certain load condition. The analysis must be 
nonlinear even if geometrical nonlinearities are not considered due to material nonlinearities asso- 
ciated with failure. These material nonlinearities enter the analysis when damage is detected and 
the local softening of the material causes the stiffnesses to change. Second, the stresses in each 
layer of the composite are calculated using a certain approach, such as the classical lamination 
theory. Then, those stresses are compared to the specific allowable strengths by way of a failure 
criterion that is used to determine the locations where failure initiates in the structure. The use of 
a failure criterion allows for the characterization of the failure; that is, it allows for a mode of fail- 
ure to be determined. Depending upon the mode of failure, a local degradation of the material 
properties is then carried out at the locations where failure is detected. There is then an iteration 
due to the need for equilibrium to be re-established. As the local stiffnesses of the structure have 
changed when the material properties are degraded, under the same pressure level, the structure 
no longer is in equilibrium. During the iteration process, generally, more failures are encountered, 
the progression of the damage can be assessed, and equilibrium re-established again. At some 
point in the iteration process, at each specific load level, no further damage is predicted, equilib- 
rium is satisfied, and the analysis proceeds to the next pressure level. 

The manner in which the progressive failure methodology was performed in [3] was as fol- 
lows. The stand-alone program used the strains computed by STAGS to calculate the stresses at 
the two interfaces of each layer (above and below the centroid of the layer). The program then 
used these stresses in the maximum stress criterion. If any of the conditions in the criterion were 
satisfied, then the mode of failure was identified; otherwise, the program returned to STAGS, the 
pressure was increased to the next load step, the strains were calculated, the program retrieved the 
new strains and, the process of analyzing failure was repeated. If conditions for failure were found 
and tensile matrix failure was identified, then the elastic modulus normal to the fibers and the 
inplane shear modulus were multiplied by a factor of p, specifically p=0.2. This level of property 
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reduction followed the work of Camanho and Matthews in [5], who followed the findings of other 
researchers, as cited in [5]. These reduced moduli were then input to STAGS, a new analysis with 
these new properties was executed to compute the strains, and the process was repeated at a spe- 
cific load level until no more failures were detected. The results obtained using this procedure 
indicated that first failure was due to matrix cracking at the clamped boundaries in all three cylin- 
ders, as expected from the results in [T,21. The first matrix failure pressures and circumferential 
locations and through thickness locations within the cylinder wall were similar to the ones pre- 
dicted by the geometrically nonlinear analysis in [1,2], The first fiber failure pressures, however, 
were lower than the pressures predicted in [1,2] using the geometrically nonlinear analysis, the 
maximum stress failure criterion, and ignoring the effects of previous matrix failures. However, 
the mode of failure was found to be due to fiber compression, as found in the previous study. 
Between the first matrix failure and the first fiber failure events, the matrix cracking pattern found 
was a function of each particular cylinder lay-up. For example, a remarkable am ount of cracking 
occurred in the axially-stiff cylinder as compared to the other two. The study concluded that no 
leakage path was formed because there was always at least one uncracked layer in each of the cyl- 
inders. 

The work of McMurray and Hyer sets the stage for the present study. They characterized the 
matrix failure scenario up to the point where the first fiber failure occurred and found that leakage 
was unlikely to develop in elliptical composite cylinders without failing some amount of fibers. 
So if the pressure is increased above the first fiber failure pressure, at some point, leakage will 
occur, and the pressure and modes of failure involved can be determined. The present study will 
continue to increase the pressure and determine when a leakage path forms and what modes of 
failure are involved when leakage occurs. However, a secondary objective of the work here will 
be to explore several issues that the failure scenarios described by McMurray and Hyer raise. 
Those issues originate from the fact that the work in [3] used a very specific progressive failure 
approach. Recall that in [3], the maximum stress failure criterion was used and that the exten- 
sional modulus perpendicular to the fibers as well as the inplane shear modulus were reduced, or 
degraded, by a factor of (5=0.2 when matrix tensile failure was detected. Legitimate questions 
arise naturally as to how the failure scenario for the elliptical cylinders would differ if different 
failure criteria and other factors of (3 and degradation schemes were to be used. Would the results 
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differ if instead of multiplying the elastic moduli by a factor of 0.2, they are multiplied by 0.05, or 
not reduced at all? After the initial failure is encountered, would degrading the moduli thereafter 
at each pressure level, thereby representing continuing softening, produce a different result than 
degrading the moduli only once? What type of impact does the degradation of the inplane shear 
modulus have on the failure scenario? That is, if the shear modulus is not degraded, but the exten- 
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criterion such as the Hashin or Tsai-Wu criteria predict more failures? If an interactive failure cri- 
terion was to predict more matrix cracks, would there still remain an uncracked layer before first 
fiber failure is detected? To answer these questions, the present work will use a similar approach 
to the one used in [3], but additional failure criteria and material degradation schemes will be 
explored and STAGS will be used to perform the progressive failure analysis. However, unlike in 
[3] that used the STAGS finite-element code to calculate strains, but utilized a stand-alone pro- 
gram to perform the progressive failure analysis, the present work takes advantage of an added 
capability to STAGS that performs the progressive failure analysis internally. Using STAGS Pro- 
gressive Failure Analysis (PFA) option, the user can choose from among five failure criteria. The 
user can also choose from a number of material degradation schemes. It is important to note that 
some of the algorithms used by STAGS are still in their experimental stage; thus, a third objective 
to this study is to scrutinize the performance of various algorithms. Specifically, at issue is the 
order of the computations in the coding within STAGS and whether or not a particular computa- 
tion can lead to unusually large numbers when moduli are degraded considerably, and whether to 
degrade engineering properties, which then degrade the stiffness matrix, or whether to directly 
degrade the stiffness matrix. 


The remainder chapters are organized so that the thesis begins by building model complexity 
and developing an understanding of the PFA in STAGS. This is done by first discussing the man- 
ner in which the simple case of a flat plate made of a single-layer of graphite-epoxy material is 
modeled and handle by STAGS using the PFA, and then by studying a single-layer circular cylin- 
der subjected to an internal radial load. Subsequently, after achieving an adequate understanding 
of the PFA in STAGS, attention turns to the pressurized elliptical cylinder problem. This is 
divided into two parts. The first part involves the comparison of matrix cracking scenarios, 
between the first matrix failure and the first fiber failure events, resulting from using various deg- 
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radation schemes and then different failure criteria for three different cylinder orthotropies; 
namely, quasi-isotropic, circumferentially-stiff, and axially-stiff. The different failure criteria to 
be compared include the maximum stress, Hashin, and Tsai-Wu criteria. In investigating different 
degradation schemes, only the maximum stress failure criterion is used. The issue of computa- 
tional algorithms is addressed at this point. The second part of this investigation assesses the pres- 
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This second part uses the maximum stress failure criterion and a particular degradation scheme. 
Finally, conclusions and some final remarks on future work are given. 
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Chapter 2 Use and Characterization of Computational 
Analysis Software 


STAGS [6] is a finite-element code mainly designed to analyze shell structures. These shells 
can be modeled in STAGS by specifying geometry, material types, loading conditions, boundary 
conditions, and finite elements that best describe the physics of a particular problem. Among the 
materials that can be modeled in STAGS are fiber-reinforced laminated composites. These com- 
posites can have up to 100 layers, with each layer having a different fiber orientation. The types of 
loading that STAGS allows shell models to be subjected to ranges from simple to complex. Sim- 
ple loadings that include prescribed displacements, point loads, line loads, surface tractions, ther- 
mal loads, and pressure loads can be combined to create specific loading requirements. The 
boundary conditions of a shell structure can be assigned in STAGS by specifying the individual 
degrees of freedom at a point, line, or surface in a shell unit to be free, fixed, or to have a pre- 
scribed value. The geometry of the structural model can be constructed by using a variety of ele- 
ments from the STAGS element library. Among the standard elements commonly used to 
construct shell structures, are planar, 4-node quadrilateral plate or curved shells. These shell ele- 
ments can have four, five, or nine integration points. Some of the standard analysis capabilities of 
STAGS include static stress, static stability, vibration, and transient analysis. These can include 
geometric and material nonlinearities. A new additional capability of STAGS is a progressive 
failure analysis, or PFA. This particular analysis, which includes geometric and material nonlin- 
earities, is used to evaluate failure initiation and progression and is the basis for analyses here. 

With the above simplified overview of some of STAGS modeling and analysis capabilities, 
the next sections will proceed to discuss in detail the specific features that will be used in this 
study, and also will validate the results of simple problems solved by STAGS. 
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2.1 STAGS Progressive Failure Analysis 

The basic steps of the PFA in STAGS can be summarized as follows: Assuming a general 
structure in which the displacements at a load step N have already been calculated and equilib- 
rium has been satisfied, begin the next load step. With an extrapolated displacement guess {A} i=0 
from the previous steps, the extrapolated strains and curvatures are calculated at the midsurface of 
the laminate, {e°} 1=0 and the {K°} i=0 , respectively. (The first time the analysis is performed, step 
N=0, a linear analysis is used to find the displacements { A} j=0 to be used as a guess for step N=l.) 

The extrapolated strains through the thickness are then computed according to the classical lami- 
nation theory, 

{ £ }; = o = {e°}, = 0 + ;{K O } < = 0 (2.1) 

where C, is the distance from the midsurface to any location within the laminate. The extrapolated 
stresses are calculated using these strains and the transformed reduced stiflness matrix from step 
N, 

{<*} = [0]{e} (2.2) 

where Q is the transformed reduced stiffness matrix of the undamaged material. A failure crite- 
rion is then evaluated using the extrapolated principal material coordinate system stresses. 
Because the extrapolated displacements {A}j = o> and consequently the extrapolated midsurface 
strains and curvatures are different, perhaps higher, than the ones in step N, the resulting extrapo- 
lated stresses may be predicted to cause failures. If this is the case, depending upon the failure cri- 
terion used and the mode of failure, the engineering properties are then degraded according to a 

particular degradation scheme, and a new stiffness matrix, [Q^] j=0 (where d stands for degraded), 
is formed. (The STAGS analysis offers two options for computing this new stiffness matrix com- 
ponents using the degraded engineering properties. However, a third option, which degrades the 
components of the stiffness matrix directly, that is, without degrading the engineering properties, 
is also possible. The details associated with these computational algorithms will be discussed in a 
later chapter and, as mentioned earlier, their performance evaluated.) After the reduced stiffness 
matrix is formed, the stresses are recomputed, that is, 

O fs),,o (2.3) 
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Also, the first variation of the energy functional, {R} i=0 , is computed using the [£^] i=0 and 
{A} /=0 . At this point a Newton-Raphson iteration is performed to determine the second variation 

or tangent stiffiiess matrix, [AT r ] i=0 , using [(/] 1=0 and {A} i=0 . With the tangent stiffiiess matrix, 

infJPmpntfll HiCTllar'pmptitQ A \ . a ran hr fminH cnhnrm riT_1. _ _ — /'EM TV* » 

r ——7 “ k—j i—v’ * U w t AV i i=U* x A1WA «* 

corrected displacement prediction can be determined by {A}^ = {A}^ + 8{A} /=0 . The strains 
and curvatures at the midsurface are recomputed using (A} {=1 and the strains through the thick- 
ness are re-calculated, becoming {e} /=] . New stresses are found using these new strains and the 
transformed stiffiiess matrix from step N, [Q]. The failure criterion is re-evaluated, the material 
properties are degraded, and the stiffness matrix is again degraded resulting in a new transformed 

stiffiiess matrix, . It is important to note that these last steps are necessary to determine the 

accuracy of the failures predicted using guessed displacements. In other words, it may be found 
that at the current iteration the displacement solution produce more, less, or no failures with 
respect to the previous iteration step. If more or less failures are found, it is due to the fact that in 
the previous iteration a lower or higher estimate of the displacements was made. Therefore, the 
engineering properties need to be degraded, or reinstated, where appropriate, and the stiffnesses 
need to be updated. If the solution is not yet converged for z=l, a new iteration is performed (z=2). 

If no failures are found, then the solution is converged, the stiffiiess matrix remains and 

this becomes [Q] (see equation Eq. (2.2)) for the next load step. 

In STAGS, the failures due to high stresses can be found by using any of five different failure 
criteria. They include the maximum stress, Hashin, Tsai-Wu, maximum strain, and Chang-Chang 
criteria. Though the first three will be used and discussed in later chapters, for a matter of conve- 
nience, these three criteria are defined here so that the reader can refer to them through out the 
remainder of this study. The failure modes for the maximum stress criterion can be written as fol- 
lows by using STAGS nomenclature and protocol, and considering a state of plane stress: 

Fiber tensile failure: 


Fiber compressive failure: 


Matrix tensile failure: 




Matrix compressive failure: 



Inplane shear failure: 



where. 


X t = tensile strength in the fiber direction 

X c — compressive strength in the fiber direction (a positive number) 
Y t = tensile strength normal to the fiber direction 


(2.5) 


( 2 . 6 ) 


(2.7) 


( 2 . 8 ) 


Y c — compressive strength normal to the fiber direction (a positive number) 

= inplane shear strength 

and, CTj, 02, and ij2 are the inplane stresses in the usual principal material coordinate system nota- 
tion. The failure modes for the Hashin criterion are denoted as follows: 

Fiber tensile failure: 



Fiber compressive failure: 



Matrix tensile failure: 



(2.9) 


( 2 . 10 ) 


( 2 . 11 ) 
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Matrix compressive failure: 



Inplane shear failure: 



The Tsai- Wu failure criterion is defined by 

F,a, +F 2 CT 2 + F 11 or^ + F 22 o| + 2F, 2 a 1 ci 2 + F 66 T2 2 > 1 
where i 


( 2 . 12 ) 


(2.13) 


(2.14) 


**1 


X c 


F 2 
F\1 


I_ 1 p _ _i_ p — _L_ 

Y t Y c 11 22 Y t Y c 

-2jF^ 2 F 66 = ± 

*xy 


(2.15) 


For the Tsai-Wu criterion the failure modes are determined as follows: 

Fiber tensile and compressive failure is assumed to occur when the following terms dominate: 


Fj<y, +Fuaf + F 12 CTiCJ2 (2.16) 

The sign of aj determines whether failure is due to tension or compression. Matrix tensile or com- 
pressive failure is assumed to occur when the following terms dominate: 

F 2 o 2 + F u a i a 2 (2.17) 

The sign of a 2 determines whether failure is due to tension or compression. The inplane shear 
failure is assumed to occur if the following term, dominates: 

Fg 6 t 12 (2-18) 


The different standard material degradation schemes that can be adopted by the PFA in 
STAGS are summarized in Table 2.1 according to failure criteria and mode of failure. 
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Table 2.1: Standard degradation schemes used by the PFA in STAGS 


Failure Criterion 

Mode of failure 

Engineering properties 
degraded 

maximum stress 

fiber tensile and 
compressive 

^l»G 12 ,vi 2 , and v 2 i, 

nr P, si nH \/~ . 

15 ■ Zi 

matrix tensile and 
compressive 

E 2> G 12 > v 12> 311(1 v 21> 
ori? 2 , v 12 , and v 21 

inplane shear 

G \2 

Hashin 

fiber tensile and 
compressive 

E b G \2’ v 12> 311(1 V 21 

matrix tensile and 
compressive 

E 2 .» G 12> v 12> 311(1 V 21 

inplane shear 

G \2 

Tsai-Wu 

fiber tensile and 
compressive 

E \> v 12 , and v 21 

matrix tensile 

1 

inplane shear 

G n 


The amount by which the engineering properties are degraded can vary. When a failure is 
detected and the failure mode is identified in STAGS, the engineering properties are multiplied by 
a factor P, where 0 < p < 1 . This factor can be either applied once after the first failure at a point 
is detected (non-recursive degradation), or the factor can be applied after first fiber failure at a 
point is detected and at each subsequent load step (recursive degradation). In other words, in a 
non-recursive degradation scheme, when a failure at a point is first detected at certain load step, 
the engineering properties are reduced by P only at that step. Thereafter, for subsequent load 
steps, the degraded properties at that point do not change. In a recursive degradation scheme, for 
subsequent load steps, the engineering properties are again multiplied by P at each load step. So 
after the third load step, for example, after failure, the engineering properties would be multiplied 

by a factor of p 3 . It is important to note that, for example, with matrix failure mode. Table 2.1 
indicates that E 2 , G 12 , v 12 , and v 2] are all degraded. The degradation of E 2 , v ]2 , and v 21 violates 
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the well-known reciprocity relation and result in an unsymmetric stiffness matrix. However, 
STAGS forces the stiffness matrix to be symmetric and E\ is not degraded. 


According to [6], the PFA in STAGS can only be performed on the standard 4-node quadrilat- 
eral plane elements from its element library. These elements, as mentioned earlier, can have four, 
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composite, because a single quadrilateral element is generally formed by several layers, the total 
number of integration points depends on the number of layers and their location within the layers. 
For a four-layer laminate, for example, assuming that a quadrilateral element with four integration 
points in a plane is chosen and that planes are located at each layer interface, the element would 
have a total of 32 integration points. At each one of these points the strains, stresses, and failure 
criterion involved in the PFA are evaluated. Before applying this analysis to the particular prob- 
lem to be investigated in this study, the following sections will examine and characterize the PFA 
in STAGS by analyzing two simple cases. 


> n rk»- n r- 




2.2 Study of Simple Cases 

To gain an insight on the actual performance of the PFA capability, a single-layer composite 
square plate under a tensile inplane load and a single-layer circular cylinder subjected to internal 
radial load will be studied. Analytical calculations will be used to determine failure loads and dis- 
placements for both cases, and these will be compared with the results computed by STAGS. 


2.2.1 Single-Layer Square Plate Under Inplane Tensile Loading 

The schematic of a perfectly flat plate can be seen in Fig. 2. 1 . In the figure only the plate’s ref- 
erence surface is depicted. The geometry of the plate is a square of sidelength a and thickness H. 
The sides of the plate are parallel with the x- and ^-directions of an orthogonal coordinate system 
located at the centroid of the plate. A uniformly distributed tensile load,/ is applied in they-direc- 

tion aty=+a/2. Also depicted in the figure are the directions of the inplane displacements, u° and 
v°, and normal displacement, w°, as well as the rotations, ru°, rv°, and rw°, as defined by the con- 


15 


ventions used in STAGS. Because of the nature of the formulation of the elements in STAGS, 
those rotations must also be specified at the boundaries. The boundary conditions at the edges of 



Figure 2.1 Geometry and nomenclature of a square plate reference surface subjected to an 

inplane tensile loading 


the plate are assumed to be as follows: 

x — u° = 0, w° = 0, ru° = 0, rv° = 0, rw° = 0, v°unspecified 

x — — 2 w ° = ru° = 0, rv° = 0, rw° = 0, n°and v°unspecified 

(2.19) 

y = + - w° = 0, ru° = 0, rv° = 0, rw° = 0, N x =/ n°and v°unspecified 
y ~ v° = 0, w° = 0, ru° = 0, rv° = 0, rw° = 0, «°unspecified 


where N x is the edge stress resultant defined by 


H 



2 


( 2 . 20 ) 


The plate is fabricated of a single-layer of graphite-epoxy with fibers oriented in the direction of 
the x-coordinate. The fiber direction is indicated in Fig. 2.1 by the three lines near the lower cor- 
ner of the plate. The engineering properties, failure properties, and the thickness of the single- 
layer are assumed to be 


16 


2Tj = 18.85 Msi 


E 2 = 1.407 Msi 


Gj 2 = 0.725 Msi Vj 2 = 0.300 h — 0.0055 inches 

X t = 200 ksi X c = 180 ksi (2.21) 

Y t = 7.25 ksi T c =180ksi 
S xy = 14.5 ksi 

These properties will be used throughout this study. The side dimension, a, for this problem is 5 
in. 


The stress-strain relations of this particular orthotropic plate for a state of plane stress can be 
written as 


( 2 . 22 ) 


where Q\\, Qi 2 > Qn* and C?66» are th® reduced stiffiiesses, which in terms of engineering proper- 
ties, are given by 

E, 


X' 


Qw Qu o 



°y 

► = 

Qn Qn o 

• < 

S' ’ 



. 0 0 0 6d 


Xxy. 


Qn = j 

Q\2 

@22 


V 12’ 1 * * * V 21 

E _2 ' v 12 
1-V, 2 -V 2 i 


(2.23) 


1 ~ v 12' v 2l 


@66 = G 12 

It is important to emphasize that due to the orthotropic nature of the material, v, 2 * v 21 . In fact, 
these ratios are related by the reciprocity relationship, 

Ej E 2 

— = — (2.24) 

V 12 V 21 

In Eq. (2.22), the only nonzero stress is the stress in the direction of the load, the y-direction. The 
stress in the fiber direction, o x , is zero because the plate is free to contract, due to the Poisson 
effect, in that direction. The shear stress, x xy> in the other hand, is zero because the uniform char- 
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acter of the applied load causes no shear deformation. Within the context of simple plate theory, 
the displacements at any point within the plate can be expressed as follows: 

u(x,y,z) = u°(x,y)-zru°(x,y) 

v(x, y, z ) = v°(x, y)-z rv°(x, y) ( 2 - 25 ) 

w(x,y,z) = w°(x,y) 

The above equation can be simplified by the facts that the plate is assumed to be perfectly flat 
before and during loading, and that the load is uniform and unidirectional. In other words, due to 
the transverse displacement, w, being zero everywhere in the plate, the only displacements to be 
consider are u and v. In addition, because the load is uniformly applied in the y-direction, the u- 
displacement is only a function of x and the v-displacement is only a function of y. As a conse- 
quence of these assumptions, the displacements in Eq. (2.25) reduce to 


u(x,y,z) = u°(x ) 
v(x,y,z) = v°(y) 
w(x,y,z) = w°(x,y) = 0 

The complete nonlinear strain-displacement relations for the reference surface are given by 

e x (x,y) = ^irf^Y + f^Y + f- 
x dx 2 l\dx) \dxJ \dx 


(2.26) 


f dw <> \ 2 “ 
di ) . 
2 . 


, s dv° ITfdu 0 ] 2 , (dv 0 } 2 (dw°\ 2 ~\ 

c y (x ’ y) = T y + Ilk* > + ^J J 

*t(x, y ) - if as'*’ +2?* 

2\_ dy dx dx dy dx dy dx 


(2.27) 



If all of the assumptions discussed above are considered, the strain-displacement relations become 

, . du° 1 f du °\ 2 

^ x) -Tx + iVTx) 


, , dv\ \(dv\ 

y (y) = Tv + ih-J 


.2 (2.28) 
, +-[— I 

dy 2 v dy 

y(x,y) = 0 

If G y from Eq. (2.22) is substituted into the total potential energy of the plate, the total potential 
energy is given by 
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n - 5 ll J*. J*. <vy*** - Mi) < 2 - 29 > 

2 2 2 

where P = f ■ a, gives the following 

H a a 

n = eps^A-Mi) (2 ' 30) 

2 2 2 


The strains e x and can be related by ct x which, as mentioned previously, is zero. Thus, from the 
first equation of Eq. (2.22), 



(2.31) 


If this strain is substituted into Eq. (2.30), and since z y is only a function of y, the total potential 
energy of the plate can be expressed as 


n = -Ha- 
2 


r Qn ' Q\\~Q \2 


Q 


11 


JV>-MD 


The first variation of Eq. (2.32) is given by 

m=Ha- 


1 Qii ■ fin ~ Qn 


Q ii 


jy V SepdM>Sv“(j) 


In terms of displacement gradients, the first variation can be written as 
5n = H-a 


' Qn ' Q\ \ - Q\2\ n 


(2.32) 


(2.33) 






The variation of + m and integration by parts of Eq. (2.34) results in 


SI1 = Ha 


f 0.22 ' 0 11 “ 012^ 
fill - 




dy 


-Pbv 


ID 


(2.35) 


P Ha- f ^ 22 ' Q 1 j- j {— +Ife°)Vl +^f)\bv°dy 

J_a ^ Q u )dy \dy 2 ^dyd J\ dy J f 
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Setting the first variation to zero yields the Euler-Lagrange equation for v°, and the corresponding 
variationally consistent boundary conditions aty=±a/2. That is. 


Ha 


Qn ' Qu ~ Q\t} 




Q 


11 


j 


d 

dy 




= 0 


(2.36) 


. . .. of a\ 

ana using me iaci mai ov -zj 


„ . of a\ . . 

u , since v 1S Known > 


Ha- 


Qji 




(2.37) 


Equation Eq. (2.36) indicates that the expression in braces is constant. From Eq. (2.37), the value 
of that constant can be determined and the equation can be written 

(£'*£TX'*S — 


2 \ 


Ha 


Qn ■ Q\\~Qn 

Qn 


H-a-E, 


(2.38) 


The above expression is a cubic equation in dv°/dy that can be solved at any load P. If the solution 
to Eq. (2.38) is denoted as dv°/d y=j{P), and this displacement gradient is integrated with respect 
toy, then the displacement v° at anyy-location for any load P is found by 

v° = /(F)y + C, (2.39) 

where Cj is a constant of integration to be solved by applying the displacement boundary condi- 
tion on v° aty=-a/2. Enforcement of Eq. (2.19), results in 

v° =AP){y + fj (2.40) 


To find the ^-displacement, Eq. (2.28) and Eq. (2.31) or Eq. (2.38) can first be used to relate the 
displacement gradients by the following equation 


du° + \fdu°\ 2 _ 0\2 f dv° + \f 2 \ 
dx 2\dx) Q u \dy 2\dy ) ) 


V 21 _P_ 1 

"" 0-g) 


(2.41) 
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If the value of dv°/dy found in Eq. (2.38) is used in Eq. (2.41), the axial displacement gradient, 
d«°/dx, can then be computed. Referring to the solution of Eq. (2.41) as du°/dx=g(/ > ), integrating 
this gradient with respect to x, and applying the displacement boundary condition on u° at x=+a/ 


2, the ^-displacement is found to be 


,,° = n(P\( r _ 

ov ~ 'V" u 




To calculate the load at which the plate will fail due to excessive tensile stress perpendicular to 
the fibers, the stress o y in Eq. (2.22) is substituted by its allowable strength, Y t , resulting in 

r.-Gu-Sa + eaV P-«) 

where the subscript t indicates failure value due to tension. Substituting Eq. (2.31) and Eq. (2.28) 
into Eq. (2.43) gives 


r ‘-— t 


(2.44) 


Recognition of the right hand side of Eq. (2.44) in Eq. (2.37), and substitution of the right hand 
side for Y t from Eq. (2.44), the failure load is found to be, from Eq. (2.37), 


+(*!)) (2.45) 

The displacement gradient at failure, (dv°/dy) t , is computed using equation Eq. (2.44). 


To complete the objective of this simple problem, which was to compare STAGS results with 
the above analytical solutions, as well as to understand the PFA, numerical results for the dis- 
placements v° and u° were obtained by modeling the square plate with a finite-element mesh com- 
posed of a total of 25 quadrilateral elements using the E410 quadrilateral shell element from 
STAGS element library. Each element was assigned with four integration points in a plane 
located at the top and bottom of the layer, for a total of eight integration points per element. A line 
load was applied at one of the edges of the mesh and boundary conditions were set so that they 
were consistent with the constraints given in Eq. (2.19). A nonlinear analysis with the PFA capa- 
bility was performed. The PFA was directed to use the maximum stress failure criterion and to 
multiply E 2 , G\ 2 , v ]2 , and v 21 non-recursively by a factor of p=0.2 for matrix tensile failure 
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mode. In the STAGS analysis, the load was increased from zero in steps of 20 lbs, except near the 
expected failure load, where the load increments were of 0. 1 lbs. The load steps were increased up 
to a load well beyond the value to cause the failure of the plate. From a physical point of view, 
there is no meaning for the results obtained after the plate failed, since for this problem every ele- 
ment had failed and the plate can not longer carry any additional load. However, for the purpose 
of understanding the PFA in STAGS, the results beyond failure, particularly the increased strains 
and displacements, were of interest as they led to a further understanding of how STAGS works. 

The plate modeled in STAGS failed at about 200.4 lbs due to matrix tension failure. The ana- 
lytical failure load predicted by Eq. (2.45) was found to be precisely that value. Figures 2.2 and 
2.3 show the comparison between analytical (solid line) and STAGS numerical (triangular sym- 
bols) displacements v° and u° at particular locations in the plate and at certain load steps. Specifi- 
cally, the figures show the v°-displacement at the edge y=+a/2 and the ^-displacement at the 
edge x=-a!2 for loads ranging from 100 to 260 lbs. Although the loads were increased from zero, 
only loads above 100 lbs are shown because more significant displacements began to take place at 
these load levels. From both figures can be seen that there is excellent agreement between the ana- 
lytical solution and STAGS results. 

In Fig. 2.2 it is shown that the v°-displacement increases with increasing load. This is 
expected since the free edge of the plate is being stretched more and more with increasing load. At 

the failure load, it can be observed that the load-v°-displacement relation exhibits a discontinuity. 
This jump is due to the sudden reduction of the extensional modulus E 2 by a factor of P=0.2 when 
failure is encountered. Note that at the failure load, E 2 in Eq. (2.38) is multiplied by 0, causing the 

displacement gradient dv°/dy to change suddenly by certain amount. As a result of that change, 
the v°-displacement in Eq. (2.39) experiences a jump at the failure load. After that jump, the mag- 
nitude of the ^-displacements are in the order of 1 0' 1 inches; while, before the failure, the magni- 
tude of the displacements are in the order of 10 inches. Also in Fig. 2.2, it can be observed that 
after failure, the slope of the load-v°-displacement curve has changed. This change in slope can be 
attributed to the fact that, before failure, the stress o y and strain e y and, therefore, the load and v°- 
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displacement, are related by the original value of the extensional modulus, E 2 . After the failure, 
however, those parameters are related by the reduced extensional modulus, p E 2 . 

Figure 2.3 shows that the ^-displacement of the edge that is free to contract due to Poisson’s 
effect increases with increasing load. This behavior is typical of a material that is being stretched 
more and more as it is pulled in one direction, causing a larger and larger contraction perpendicu- 
lar to the direction of the pull. The fact that the M°-displacements in the figure are positive indicate 
that the free edge moves in the positive direction toward the fixed edge. The small discontinuity in 

the load-iAdisplacement curve which can be observed in the figure is associated with the change 
in the displacement gradient dv°/dy mentioned above. Note that since immediately after failure E 2 
is multiplied by P and v 21 is also multiplied by p, to maintain the reciprocity relation of Eq. 

(2.24), the term dv°/dy in the denominator in the right hand side of Eq. (2.41) is the only term 
responsible for the sudden change in d«°/dx, which in turn produces the discontinuity in the load- 
tAdisplacement curve at the failure load. The magnitude of the ^-displacements are in the order 
of 1C 4 inches. 

From this simple case of a single-layer square plate under inplane tensile loading, it was deter- 
mined that the numerical results given by STAGS and its PFA capability are in good agreement 
with the analytical solutions described. The errors in the results, which were found to be insignif- 
icant, are believed to be due to the computational algorithm, which is generally expected of finite- 
element codes. In the next section, a slightly more complicated case with different geometry and 
loading will be study analytically and using STAGS. 


2.2.2 Internal Radial Loading of a Single-Layer Circular Cylinder 

Seen in Fig. 2.4 is a diagram describing the next problem designed to study the PFA option in 
STAGS. The geometry of the reference surface of a single-layer circular cylinder of radius R and 
length L subjected to a uniformly distributed radial load, denoted by p Q , is depicted. A global ref- 
erence system is located at the geometric center of the cylinder. The coordinate C, is measured 
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Figure 2.4 Geometry and nomenclature of a circular cylinder reference surface subjected 

to an internal radial loading 

from the middle of the single-layer, like z in Fig. 2.1. The axial, tangential, and outward displace- 
ment directions are indicated by u°, v°, and w°, respectively. The boundary conditions are given 

u° = 0, ru° = 0, rv° — 0, v°, w°, and rw° unspecified 

(2A6) 

ru = 0, rv = 0, u , v , w°, and rw° imspecified 

Because of the formulation of the element in STAGS, the following other conditions were 
imposed: 



u°, w°, rv°and rw ° unspecified everywhere else 
ru° = 0 everywhere 

v° = 0 at specific locations to prevent rigid body motion 


(2.47) 


The cylinder is made of a single-layer of graphite-epoxy, with the same material properties 
and thickness as in the previous case. A list of those was given in Eq. (2.21). The radius of the cyl- 
inder is 5 inches and the length is 12.5 inches. 
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If a state of plane stress is assumed, as it was done for the case of the plate, Eq. (2.22) through 
Eq. (2.24) remain valid for this case. Within the context of simple cylinder theoiy, the three com- 
ponents of the displacement, analogous with Eq. (2.25), can be written as 


0,0 = u°(x, 0 ) - C, ru°(x, 0 ) 

v(x, 0, 0 = v°(x, 0) - ; rv°(x, 0) (2-48) 

w(x, 0 , Q = w°(x, 0 ) 

With boundary conditions, Eq. (2.46), this circular cylinder under the uniform loading p 0 is free to 
expand in the outward direction, is free to contract axially, and the response of the cylinder is 
expected to be axisymmetric. Therefore, the following assumptions can be made. 

v(x, 0, Q = 0 
— ( ) = 0 

8Q y } (2.49) 


— ( ) = — ( ) 
dx y ’ dx K ' 

Simplifications of Eq. (2.48) can be made by using Eq. (2.49). As a result, Eq. (2.48) reduces to 


u(x, 0, Q = u°(x ) 

v(x,0,Q = 0 (2.50) 

w(x, 0, Q = w°( a constant) 

The complete nonlinear strain-displacement relations for the reference surface in cylindrical coor- 
dinates are derived from [9] 


^-^tR[© 2+ (gy + (£T] 

^■ 6 > ■ + *) 2+ G¥-$ 2 ] 

v , m _ \\\du° ,du°(\du°\ dv° dv°(ldv\w°) dw°(ldw° v°)1 

Y ^’ 0) ' i\m + Tx^mJ + d-x + T) + di \im ~r)\ 


(2.51) 
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If Eq. (2.50) is substituted into Eq. (2.51), the reference surface strains simplify to 

o (\2 


s x (x, 0) = 


If**)' 

2 \dxJ 


du 

dx 2 Wjc 
0 




(2.52) 


7*e(*» 0) = 0 

Since the only nonzero stress is the circumferential stress, ct 9 , the total potential energy for the 

t 

radially loaded circular cylinder is given as [10] 

H L L 

n = \\\ jj* |\ a^dxRMdC, - x J>j(vv”y? + 1 w ° l ) dx (2.53) 

2 2 2 

The circumferential stress, a e , in Eq. (2.53) can be written in terms of strains using the second 
equation in Eq. (2.22), with the appropriate coordinate change. This results in 
Hi L 

n = 71 J 2 // ?l (£12 ' e x + 022 • 8 e> 8 e • RdxdC, - n J 2 £ 2p 0 (w°R + lw° 2 ^ dx (2.54) 
~2 ~2 ~2 

Since a x is zero, the cylinder is free to contract axially as it expands outward. Thus z x is related to 
the circumferential stress e 9 by 

012 

8 * = ’ 8 e (2-55) 

Substitution of Eq. (2.55) into Eq. (2.54) and use of the fact that the circumferential strain is not a 
function of allow Eq. (2.54) to be rewritten as 
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The first variation of the total potential energy is given by 
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(2.57) 


If the circumferential strain in Eq. (2.57) is expressed in terms of the outward displacement, as in 
Eq. (2.52), and the second term is rearranged, the first variation becomes 
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If the first variation is set to zero, the following Euler-Lagrange equation is found: 
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(2.59) 


Equation (2.59) is an algebraic equation which allows for the determination of the outward dis- 
placement as a function of the internal load. To find the axial displacement, u°, as a function of the 
axial coordinate and the internal pressure, Eq. (2.55) expressed in terms of displacements can be 
used. That is. 


du 

dx 


— »-(t+Kt)’) < 2 -«» 


The Euler-Lagrange equation found in Eq. (2.59) can be used to rewrite the above equation as 

2 ~Po R Q\2 _ V 21 Po R 
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H 


(2.61) 


From this equation, the displacement gradient, dw°/dx, can be computed and used to determine the 
axial displacement, u°. If the solution to Eq. (2.61) is denoted as du°/dx=h(p 0 ), this displacement 

gradient integrated with respect x, and u° is set to be zero at x=-L/ 2, according to Eq. (2.46), the 
axial displacement at any axial location and at any load level is given by 

u° = A<p 0 )(x + D (2-62) 

The load at which the cylinder is expected to fail occurs when the stress in the circumferential 
direction reaches the allowable stress. Mathematically, this is equivalent to writing 

Y t = Q\ 2 ' z xt + Q 22 ' Ht (2.63) 

where again the subscript t in the strains indicates failure values due to tension. Equation (2.55) 
and Eq. (2.52) can be used to write the axial strain in terms of the circumferential strain and the 
circumferential strain in terms of outward displacement, respectively. Thus, the above equation 
yields 
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(2.64) 




Substitution of the second factor in Eq. (2.64), which was evaluated by the Euler-Lagrange equa- 
tion, Eq. (2.59), allows for the solution of the failure load in terms of the allowable strength. That 
is, 

Y t -H 


Pot = 


R 


(2.65) 


STAGS numerical results for the failure load and the outward and axial displacements for this 
case were obtained by constructing a finite-element mesh qf the circular cylinder composed of 10 
elements in the circumferential direction by 30 elements in the axial direction. The type of ele- 
ments used was the E410 quadrilateral shell element with four integration points at the top and 
bottom of the layer, for a total of eight integration points. An internal radial loading was applied 
and boundary conditions were set so that no rigid body motion was allowed and the cylinder was 
free to expand outwardly and one end of the cylinder was free to contract axially. The load was 
increased from zero up to 13 psi in steps of 2 psi, except near the expected failure load, where the 
load was increased in steps of 0.1 psi or less. As in the case of the plate problem, the load was 
increased above the failure load only to characterize STAGS PFA capability, since, physically, 
this cylinder would lose its structural integrity as soon as the failure load is reached. The settings 
used in the PFA were the same as in the previous plate case. That is, failure was determined using 
the maximum stress failure criterion and the degradation scheme used reduced non-recursively 
the elastic moduli E 2 , G\ 2 , v i2> and v 2 ] by (3=0.2 for matrix tensile failure mode. 


The axial, circumferential, and outward displacements were recorded for each node at each 
load step. Figures 2.5 and 2.6 show the load vs. the outward displacement, w°, for any axial loca- 
tion (red triangular symbols), and the axial displacement, u°, at the free end of the cylinder (blue 
triangular symbols), respectively. The circumferential displacement was not plotted since, as 
expected, due to the axisymmetry of the problem, it was zero. Shown in the figures are also the 
analytical solutions for the displacements at each load level (solid line). As it can be observed 
from the figures, analytical solutions and STAGS results are in close agreement. 
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Figure 2.5 shows that the outward displacement, w°, increases as the internal load is 
increased. Also shown in Fig. 2.5 is that at the failure load, a discontinuity in the outward dis- 
placement occurs. The change in displacement is directly associated with the fact that the degra- 
dation experienced by the extensional modulus, P E 2 , at the failure load produces an abrupt change 

in the in circumferential strain. This sudden rise in strain causes the displacement w° to also 
increase abruptly by the second strain-displacement relation in Eq. (2.52). The magnitude of the 

normal displacements after failure are in the order of 10* 1 inches, while the magnitude of these 

displacements before failure are in the order of 10' 2 inches. Another observation that can be made 
from the load-displacement curve in Fig. 2.5 is the change in its rate of increase before and after 
the failure load. The cause of this change can be attributed directly to the degradation of E 2 - 
Before the failure load is reached, the extensional modulus in the load-displacement relation 

given in Eq. (2.59) has its original value. This results in particular values of ^-displacements that 
increase at a certain rate with the load. Immediately after the failure load is reached, E 2 is reduced 
by the factor p. As a result, the outward displacements increases at a lower rate with load. 

In Fig. 2.6, it is seen that the magnitude of the ^-displacement increases with increasing load. 
The negative signs indicate that at a given load the cylinder’s free end moved toward the other 
end, which is fixed axially. The increased shortening of the cylinder is the result of the Poisson’s 
effect becoming larger, as higher loads force the cylinder to expand more in the outward direc- 
tion. This situation is analogous to the one in the previous case in which the plate experienced a 
larger and larger contraction in one direction as it was pulled more and more in the transversal 

direction. One interesting observation that can be made in Fig. 2.6 is that unlike the ^-displace- 
ment in the plate problem, the axial displacement is continuous. To see why this is the case, Eq. 
(2.61) should be carefully examined. Immediately after failure, E 2 and v 21 are multiplied by p. 
This causes the P’s to cancel from Eq. (2.61). Therefore, the right hand side of the equation 
remains constant, and consequently the solution to the equation du°ldx~h{p 0 ) also remains con- 
stant 
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Pressure, p 0 psi 







2.3 Summary 

In this chapter the basic features of the finite-element code STAGS which will be used 
throughout the remainder of this study have been introduced. It was indicated that this code is 
generally used to perform stress analysis in thin shells. The specific feature, the PFA option, in 
STAGS was particularly discussed in detail. In the discussion, it was pointed out that the PFA 
capability allows an analysis of failure initiation and progression in a given structure. This type of 
problem forms the basis of the study to be presented in this work. 
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Chapter 3 Matrix Failure Scenarios in Internally 
Pressurized Elliptical Cylinders 
Using Different Degradation Schemes 


In this chapter the effects of different degradation schefnes on the matrix failure scenarios for 
an internally pressurized, elliptical, composite cylinder with different material orthotropy will be 
investigated. As mentioned in the introduction, a secondary objective of this study is to answer 
some of the questions raised by the work of McMurray and Hyer [3]. Recall that in their study 
they used a progressive failure analysis to study failure of an internally pressurized, elliptical, 
composite cylinder with quasi-isotropic, circumferentially-stiff, and axially-stiff material 
orthotropies. In their progressive failure methodology, they used the maximum stress failure crite- 
rion and a degradation scheme that reduced the extensional modulus perpendicular to the fibers 
and the inplane shear modulus by 80% when matrix tensile failure was encountered. Coinci- 
dently, up to the pressure levels they considered, only matrix tensile failures were encountered. 
The fact that McMurray and Hyer used such a specific failure criterion and degradation scheme 
prompted several questions in regard to how their failure predictions would change if other crite- 
ria and degradation schemes were used. Specifically, the chapter will address the following ques- 
tions: What is the effect of degrading the inplane shear modulus? How much does the factor by 
which the engineering properties are reduced influence the failure scenarios? Does it make any 
difference to reduce the engineering properties once when failure is first detected rather than to 
recursively reduce the properties at each subsequent pressure level? To study the influence of the 
different degradation schemes, the PFA option will be configured to use the maximum stress fail- 
ure criterion. 

By way of a brief outline, the chapter will be organized to begin by describing the geometry, 
nomenclature, material specifications, and finite-element model of the elliptical cylinder to be 
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investigated. Next, a method of illustration to visualize the failure scenarios will be presented. 
The chapter will then discuss the failure scenario at several pressure levels for a cylinder with 
quasi-isotropic material orthotropy resulting from using a degradation scheme that degrades non- 
recursively the extensional modulus in the direction perpendicular to the fibers, E 2 , the inplane 
shear modulus, G 12 , the major Poisson’s ratio, v ]2 , and the minor Poisson’s ratio, v 21 , when 
matrix tensile failure is detected, (see Table 2.1) This degradation scheme will be referred to as 
the baseline scheme. The results obtained for the quasi-isotropic cylinder will then be used to 
compare with the resulting matrix failure scenarios for cylinders with circumferentially-stifF and 
axially-stiff material orthotropy; the same baseline degradation scheme was used for all three 
material orthotropies. In keeping with the procedure in Ch. 2, a value of fJ=0.2 will be used. From 
that comparison some conclusions on the effects of material orthotropy on the failure behavior of 
elliptical cylinders will also be made. Finally, other degradation schemes will be used, including 
using other values of P, and comparisons between the different schemes will be made. 


3.1 Problem Formulation 

To provide a basis to study and compare the results obtained in this study, an internally pres- 
surized elliptical cylinder model similar to the one used by McMurray and Hyer is adopted. It is 
believed that the results for an elliptical cylinder can be extended to more general non-circular 
cylinders due to their common trait, namely, the varying curvature of their cross sections. Figure 
3.1 describes the elliptical cylinder by depicting the geometry of its reference surface and the 
nomenclature to be used. As can be seen from the figure, the length of the cylinder is represented 
by L and is oriented in the ^-direction of an orthogonal reference system located at the geometric 
center of the cylinder. The major and minor diameters of the elliptical cross section are designated 
by 2 a and 2b, respectively. The ellipticity of the cylinder is denoted by e and is defined here as b/ 
a. The parameter s represents the circumferential coordinate measured counterclockwise from the 
top, or crown, of the cylinder, and C denotes the circumference of the cylinder cross section. The 
wall thickness is denoted by H and the internal pressure is indicated by p a . The crown, keel, and 
sides of the cylinder are identified by the locations s/C- 0, s/C-0.5, and s/C=± 0.25, respectively. 
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Though they will not be discussed explicitly, also shown in fig 3.1 are the directions of the 


displacements and the rotations of the reference surface. The displacements are denoted by u°, v°, 
and vv°, and represent the axial, circumferential and outward normal displacements, respectively. 


The corresponding rotations about those directions are denoted as ru°, rv°, and rw°. The displace- 
ments and rotations are all functions of the jc and s coordinates. The boundary at one end of the 
cylinder is assumed to be completely clamped, while the other end is clamped but allowed to 
move freely in the axial direction, as if it was fastened to a bulkhead. Mathematically, the bound- 
ary conditions are expressed as follows: 


J i O . O n O /~\1 ® g\ 

x = -- u = 0, v = 0, w = 0, rv = 0 
jc = u = A, v = 0, w = 0, rv = 0 


(3.1) 


where A is given by enforcing equilibrium in the axial direction; that is 

fcNjb-x-a-b-p, (3.2) 

where N x is the axial force resultant defined as 

H 

K = (3.3) 

2 

The cylinder is modeled assuming graphite-epoxy material with properties as listed in Eq. (2.21). 
Three different balanced symmetric laminates are considered; namely, quasi-isotropic 
[±45 /90 / 0] s , circumferentially-stiff [±45 / 90 2 / 0 1/2 ] J , and axially-stiff [±45 / 0 2 / 90 1/2 ] j 
laminates. The fiber orientations for all the layers in these laminates are measured with respect to 
the axial direction. It should be mentioned that each laminate has at least one layer with fibers in 

the 0° direction, one layer with fibers in the 90° direction, and all have ±45° layers on the inside 
and outside. In this study laboratory scale elliptical cylinders were studied so that, eventually, 
experiments can be performed to compare with the results obtained here. From a practical point of 
view, small scale cylinders are easier to fabricate by hand using elliptical mandrels than large 
scale cylinders, and they are less expensive to fabricate and test. The dimensions of the cylinder 
are as follows: 

L = 12.5 in. a - 5 in. b = 3.5 in. 
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Figure 3.1 Geometry and nomenclature of elliptical cylinder 
3.2 Finite-Element Model 


The finite-element mesh of the elliptical cylinder is shown in Fig. 3.2. This mesh is con- 



Figure 3.2 Finite-element mesh of elliptical cylinder 

structed using five shell units with a total of 3500 elements, there being 70 elements in the axial 
direction and 50 elements in the circumferential direction. The shell elements are the plane-stress 
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E410 quadrilateral element from the STAGS element library. The elements are configured to use 
four integration points at the top and bottom of each layer of composite material. That is, for the 
eight and nine layer composite laminates studied, there are 64 or 72 integration points per ele- 
ment, respectively. The mesh spacing is uniform in the circumferential direction, however the 
mesh density in the axial direction, as Fig. 3.2 shows, increases near the clamped ends to accu- 

rnfplv HpcrriVip tlipsp *7nr»pR where hi ah axial oraHients are evnerted 

J — - ^ — 0 — — r 


3.3 Method of Illustration 


To illustrate the resulting failure scenarios here and throughout this study, the method of illus- 
tration shown in Fig. 3.3 will be used. A quadrant of the cylinder is isolated from the rest of the 



Figure 3.3 Method of illustrating one quadrant of the elliptical composite cylinder 


cylinder, as highlighted in the top left portion of the figure. This quadrant is then flattened and 
expanded in the thickness direction, which is greatly exaggerated. These two steps are identified 
in the figure. A three-dimensional map is then created, as shown in the bottom left portion of Fig. 
3.3. In this map, the non-dimensional range of the quadrant selected extends from 0 <*/£,< +0.5 
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in the axial direction and —0.25 < s / C < +0.25 in the circumferential direction. The non-dimen- 
sional range in the thickness direction is —0.5 < C, /H < +0.5 . Where Q=Q at the geometric midsur- 
face of the cylinder wall. The layers are numbered so that the inner, or bottom layer, is layer 1 . 
The finite-elements consisting of all the layers fill this three-dimensional space. Although failures 
can occur at each integration point in a layer in an element, when a matrix tensile failure is 
encountered, for examnle. a sinele red circular snhere is nlotted at th<» rantmiri of the laver in the 

* ' A A ' ■ - ■ * V " 

element in which failure of any of its integration points occurred. This procedure is done due to 
the limited resolution of this graphical representation. As the pressure is increased and more 
matrix failures are detected, more spheres are added to the map, exhibiting a clearer matrix crack- 
ing pattern. It is important to emphasize that due to the presence of the bending stiffnesses Dj 6 
and D 26 , there is not quarter symmetry in the response. However, using symmetric and anti-sym- 
metric arguments, the cracking pattern for the whole cylinder can be envisioned, and displaying 
one quadrant of the cylinder as in Fig. 3.3 results in a less cluttered display of failure locations. 
The next section will present and discuss in detail the matrix failure maps at specific pressure lev- 
els for the quasi-isotropic cylinder resulting from the use of the baseline degradation scheme. 


3.4 Failure Initiation and Progression in a Quasi-Isotropic Elliptical Cylinder 
Using the Maximum Stress Failure Criterion and the Baseline 
Degradation Scheme with P=0.2 

For the quasi-isotropic cylinder, the pressure was increased in increments of 10 psi or less to 
determine initial failure and to continue monitoring for failures until first fiber failure occurred. 
Before performing the analysis, it was determined that below 100 psi the elliptical cylinder was 
unlikely to fail. That was a rough estimate based on a circular cylinder with a radius equal to the 
average of the semi-minor and semi-major diameters of the elliptical cylinder. The maximum 
stress failure criterion was selected in STAGS to determine failure. The conditions for failure 
using this criterion were listed in Eq. (2.4) through Eq. (2.8). The PFA option in STAGS was set 
to use the baseline degradation scheme using a factor of p=0.2 when matrix tensile failure was 
detected. Matrix failures due to shear in the x-s plane could be accommodated, but up to first fiber 
failure, surprisingly, none were encountered. The output given by STAGS listed the mode of fail- 
ure, element number, layer number, and integration points where failures were detected. With that 
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information, the three-dimensional map described in the previous section was filled to display the 
failure pattern in the cylinder. 

Figure 3.4 shows the failure maps for the quasi-isotropic cylinder at four different pressure 
levels. The pressure values are indicated at the upper left hand comer of each map. The fiber 
angle 2 nd layer number are written next to each layer The red spheres in all of the failure maps 
indicate that the mode of failure at these pressures is matrix tension, or matrix cracking, which 
means excessive a 2 stress in tension in the direction perpendicular to the fibers. At the top left 
comer of the figure the failure map corresponding to a pressure of 140 psi is shown. This pressure 
is slightly above the pressure to cause first matrix tensile failure. The pressure at which the first 
failure occurred was found to be approximately 138 psi. As can be seen from the map, at 140 psi 

two elements in the inner +45° layer, or layer 1, exhibit matrix tensile cracks. The axial location 
of these failures coincides with the end of the cylinder. The circumferential locations of these 
cracks are partway between the center of the crown and one of the sides, at s/C«-0.13. It should 
be pointed out that there are no companion failures at s/C*+0.13 due to the nonzero stiffnesses 
D I6 and D 26 . These stiffnesses couple bending and twisting, producing an anti-symmetric failure 
pattern. In other words, D ]6 and D 26 cause companion failures to the ones shown at s/C~+0.37 in 
the lower quadrant at the x/L=+ 0.5 end of the cylinder. At the other end of the cylinder, at x/L— 
0.5, the matrix cracks would appear at s/C«+0.13 and at s/C«- 0.37. The first matrix failure pres- 
sure as well as the location of the first matrix cracks agreed with the failure pressure and location 
found by McMurray and Hyer in [3], where a stand-alone pre/post processing program to perform 
the progressive failure analysis was used in combination with STAGS. This agreement was 
expected because, in their work, they used a similar finite-element mesh and the progressive fail- 
ure analysis is not necessary up to the point of initial matrix cracking. The initial matrix cracks 
shown in the map corresponding to the 140 psi pressure level are believed to occur at those loca- 
tions mainly due to large values of shear stress resultant and high bending effects that originate 
from the clamped boundary conditions which forces any circumferential and normal displace- 
ments to be zero at the boundary. These suspicions agree with the findings in [1,2]. Recall, as 
stated in the introduction, that those studies showed that the maximum shear force resultants 
occur near the clamped boundaries, part way between the ends of the minor of the major diame- 


39 


ters in the four quadrants at each clamped end of the cylinder, and that the maximum bending 
moment resultants occur at the clamped ends. 


At 170 psi, shown at the top right comer of Fig. 3.4, it can be seen that more matrix tensile 
failures occur in the inner ±45° layers, and new failures occur in the outer 90° and +45° layers. 


compared to me ihu psi map, uic in uic ' *+ j icty^i udv& ui&m iy |n> > csscd circ umf breii— 

tially toward the crown and the side. Although difficult to notice, some axial progression toward 
the center of the cylinder also occurs in this layer. This becomes more clear at 200 psi by noticing 
what seems like the elongation of the round symbols. Due to the perspective of this graphical rep- 
resentation, as cracking moves axially and a symbol is added, it overlaps the symbol behind it, 
giving the appearance of an elongated symbol. Also at 200 psi, it can be observed that all cracking 
patterns continue to expand circumferentially. New cracks also have developed at the sides in the 


other 90° and -45° layers. The two 0° layers do not exhibit any cracks. At 230 psi, it seems that 
more axial progression of the cracks has occurred, in particular, near the sides. Circumferential 


progression of the cracks in the inner ±45° layers, as seen by the overlapping near s/C=0.05, begin 
to create a contiguous path between those two layers. However, because the remainder layers do 
not possess any cracks at that axial and circumferential location, there is no leakage path. More- 
over, it was determined that up to the pressure to cause first fiber failure (not shown), no complete 

leakage path had developed because the two 0° layers remained free of cracks. This findin g coin- 
cides with McMurray and Hyer’s prediction in [3]. The first fiber failure pressure was found to be 
approximately 246 psi. The amount of damage in the cylinder at that pressure was assessed by the 
ratio of the number of integration points with damage (matrix and fiber) and the total number of 
integration points. That ratio indicated that 5.5% of the cylinder had damage at the first fiber fail- 
ure pressure. The first fiber failure mode was found to be fiber compression and it occurred in the 


outer +45° layer near s/C=-0.13. This will be discussed later. One interesting point that should be 
mentioned is that the overall cracking pattern for the quasi-isotropic cylinder using the PFA 
option in STAGS was found to be similar to the failure scenario found by McMurray and Hyer. 
However, the first fiber failure pressure was found to be slightly higher in [3] than the one found 
here. It is believed that the reason for the discrepancy is that in [3], a different convergence crite- 
ria was used in the iteration process of the progressive failure analysis. 
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Figure 3.4 Matrix cracking scenario for the quasi-isotropic [±45/0/90] s elliptical cylinder at 
different pressure levels, maximum stress criterion, baseline degradation scheme, {3=0.2 


The same analysis performed on the quasi-isotropic elliptical cylinder was repeated for the 
circumferentially-stiff and axially-stiff cylinders. The matrix cracking maps for these cylinders 
can be seen in Fig. 3.5. The first column in the figure shows the failure maps at 140, 170, 200, and 
230 psi previously shown in Fig. 3.4 for the quasi-isotropic cylinder. The second and third col- 
umns show the cracking patterns, at the same pressure levels, for the circumferentially-stiff and 
axially-stiff cylinders, respectively. Figure 3.5, clearly shows how material orthotropy influences 
the matrix cracking patterns in an elliptical cylinder. It should be mentioned again that the circum- 
ferentially-stifif and axially-stiff cylinders are made of 9-layers with [±45 / 0 2 / 90, /2 ] j and 

[±45 / 90 2 / 0 )/2 ] 5 fiber orientations, respectively. Because of space limitations, the layer num- 
bers and fiber angles are omitted in the subfigures. 
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The middle column of Fig. 3.5 shows the cracking maps for the circumferentially-stiff cylin- 
der. At 140 psi can be seen that matrix cracks occur at the sides of the cylinder rather than partway 
between the crown and the side, as is the case for the quasi-isotropic cylinder. At that pressure 

level, the circumferentially-stiff cylinder has cracks in the outer ±45° layers, layers 8 and 9, the 

two outer 90° layers, layer 6 and 7, as well as the 90° layer just inside the midsurface, layer 4. 
This cylinder exhibits many more cracks than the quasi-isotropic cylmder at 140 psi because the 
first matrix tensile failure in the circumferentially-stiff cylinder occurred at a pressure of approxi- 
mately 105 psi. So a number of cracks have developed between 105 and 140 psi. The first matrix 
failure pressure and location for this cylinder, as expected^ also agreed with the failure pressure 
and location found in [3]. It is believed that the cracks initially occur at the sides of the circumfer- 
entially-stiff cylinder because of high values of the circumferential bending moment and force 

resultants there, and low resistance to cracking in the axial direction of the 90° layers. There is 

only one 0° layer to directly resist axial loads. The map associated with the pressure of 170 psi 
shows that the cracking has mainly progressed axially, but also there is some progression circum- 
ferentially. This map for 170 psi also shows that several more matrix failures occurred in the inner 
three layers. 

At 200 psi there is even more axial and circumferential progression. However, no cracks are 
seen in the 0° layer, the central layer, layer 5. At 230 psi, the progressive pattern continues, now 

causing some overlap of matrix cracks in the inner ±45° layers, layers 1 and 2, near s/C=+0.05, as 
was the case for the quasi-isotropic cylinder. In addition, new cracks developed at the side s/C— 

0.25 in the inner -45° layer, clearly forming a potential continuous path for the pressure to escape. 

However because the single 0° layer remains crack-free, the formation of a leakage path is 
impeded. Interestingly, up to the pressure to cause first fiber failure, which was found to be 253 

psi, no complete leakage path forms. It is believed that the single 0° layer remains uncracked 

because the 90° layers protect that layer. The percent damage in the cylinder when first fiber fail- 
ure occurred was approximately 6.9%. The mode of the first fiber failure was found to be fiber 

compression and the location of this fiber failure was in the outer +45° layer. More will be said of 
the first fiber failure later. 
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The furthest right column of Fig. 3.5 shows the maps for the axially-stiff cylinder. At the pres- 
sure levels of 140 and 170 psi, it can be seen that the axially-stiff cylinder exhibits a very similar 
cracking pattern to the one found for the quasi-isotropic cylinder. Even the first matrix tensile fail- 
ure was found to occur at practically the same pressure as for the quasi-isotropic case; precisely, 
the first matrix tensile failure pressure occurred at about 139 psi in the axially-stiff cylinder. Once 
again, the failure pressure and locaticn agreed with the findings in [3]. It is suspected that the fail- 
ures, as in the quasi-isotropic case, occur due to high bending effects at the clamped ends of the 
cylinder and high values of the shear stress resultant produced by the boundary condition on 
restricting the normal and circumferential displacements. At 200 psi, the cracking pattern changes 
significantly with respect to the quasi-isotropic case. A large amount of cracking occurs within 
the range from 170 to 200 psi. In fact, because of a large number of cracks that developed 
between 180 and 190 psi step (recall that the analysis was performed in 10 psi steps), the number 
of iterations allowed in STAGS had to be increased. Recall that the PFA option at a given load 
step performs a Newton-Raphson iteration to determine if more, less, or no failures occurred with 
respect to the previous iteration. If more failures occur at certain iteration, then the engineering 
properties are degraded, stresses computed, failure checked, and if more failures are detected, 
another iteration takes place. This is repeated until no more failures are found and the analysis is 
said to have converged for that load step. STAGS allows a maximum of 22 of these iterations at a 
given load step. If this number is exceeded, the program cuts the load step and repeats the analysis 
at a lower load step. The number of cuts is set by the user. In the analysis performed here, that 
number was set to 5. For the axially-stiff cylinder, the maximum number of iterations needed for 
the analysis to converge at the load step between 180 and 190 psi was exceeded and the program 
cut the load step once in half. The analysis converged at the load step between 180 and 185 psi. 
However, when the pressure was increased from 185 to 190 psi, the analysis did not converged 
even after 5 cuts in the pressure step. A special option in STAGS which extends the maximum 
number of iterations allowed to 40 was then executed. With this option, no cuts in the 10 psi pres- 
sure steps were needed. 

Also in the cracking map at 200 psi in Fig. 3.5, it can be seen that most of the matrix cracks 

are concentrated in the 0° layers, layers 3, 4, 7, and 8, near the crown, away from the clamped end. 
It can also be observed that at 200 psi the cracking has progressed quickly axially towards the 
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midspan of the cylinder. It is believed that this cracking pattern develops due to large values of the 
circumferential force resultants. Because the maximum circumferential force resultants occur 
partway between the clamped end and the midspan of the cylinder in the crown region, as indi- 
cated by the previous studies [1,2], and because there is only one 90° layer in this cylinder wall to 

resist those force resultants, layer 5, cracking in the 0° layers quickly progresses towards the mid- 
span of the cylinder in the crown area. 

j 

To actually see in the cracking map that there is an axial distance between the end of the cyl- 
inder and the circumferential line of cracks closest to x/L=f 0.5 in layer 7, a 0° layer, for instance, 
a line parallel to the s/C-axis can first be drawn through the centroid of that layer at x/L=+ 0.5. A 
second line also parallel to the s/C - axis can then be constructed that passes through the centroids 
of the round symbols forming this line of cracks. Note that the centroids of the round symbols in 
this line have to be estimated, since other rows of symbols overlap the first row. Finally, a third 
line, this time parallel to the c-axis and long enough to intersect the previous two, can be traced 
through the centroid of layer 7 at s/C=+0.25. The fact that there is a finite distance between the 
intersection point of the first with the third line and the intersection point of the second with the 
third line, indicate that the first row, and therefore subsequent rows of cracks, occur away from 
the end of the cylinder. This procedure can be followed to estimate the distance between the end 
of the cylinder and any axial location. Also at 200 psi it is shown that the only uncracked layer is 
layer 5, or the 90° layer. 

The map at 230 psi shows an impressive amount of damage. Cracks have progressed axially 
as well as circumferentially and encompass a significant percentage of the three-dimensional 

space. The layers that exhibit the most cracks are the 0° layers. As mentioned above, because 
there is a small number of fibers in the circumferential direction to resist the large circumferential 

force resultants, these 0° layers become the most vulnerable to cracking. Although difficult to see 

from the map, the 90° layer, surprisingly, remains crack-free regardless the high number of matrix 
tensile failures in the other layers. At the first fiber failure pressure (not shown), which occurs at 

approximately 251 psi, this 90° layer does not exhibit any cracks, even though 26.5% of the entire 
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cylinder was found to be damaged. Because up to first fiber failure, the 90° layer remains 
uncracked, the axially-stiff cylinder does not have a leakage path for the pressure to escape. 
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Figure 3.5 Matrix cracking scenarios for elliptical cylinders with different material orthot- 
ropy, maximum stress criterion, baseline degradation scheme, p=0.2 
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Tables 3.1 and 3.2 summarizes the pressures for first matrix tensile and first fiber failures for 
each of the cylinders. In Table 3.2, the percent damage associated with the first fiber failure pres- 
sures are also listed. 

Table 3.1 : First matrix failure pressures for maximum stress criterion 


Cylinder 

First matrix failure 


nroccnrp noi 

V/A. U1U U 


quasi-isotropic 

138 

circumferentially-stifF 

105 

axially-stiff 

139 


Table 3.2: First fiber failure pressures for maximum stress criterion, 
baseline degradation scheme, (3=0.2 


Cylinder 

First fiber failure 

Percent 

orthotropy 

pressure, psi 

damage 3 

quasi-isotropic 

246 

5.5 

circumferentially-stifif 

253 

6.9 

axially-stiff 

251 

26.5 


a. Number of failed integration points / total number of integration points 


The key results for the three cylinders just discussed are summarized as: 

i. All of the cylinders were found to exhibit the same type of failures at the indicated pressures. 
Specifically, matrix tension, or matrix cracking, which is associated with excessive stress in the 
direction perpendicular to the fibers, occurred for all pressure levels up to first fiber failure pres- 
sure. 

ii. The first matrix cracks occurred, for all intents and purposes, at the same pressures and loca- 
tions for the quasi-isotropic and axially-stiff cylinders. That pressure was close to 140 psi. The 
matrix cracks at that pressure were found, mainly, to be due to high shear stress resultants and 
bending effects. For the circumferentially-stifF cylinder, the first matrix failure occurred at a 
lower pressure than for the other two, i.e., at approximately 105 psi. The cracks in this cylinder 
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occurred at the sides and were mainly due to large values of the axial bending moment and force 
resultants. 

iii. The quasi-isotropic cylinder sustained the least amount of damage of all three cylinders, i.e., 
5.5% damage at first fiber failure pressure. The circumferentially-stiff cylinder exhibited slightly 
more cracks than the quasi-isotropic cylinder. The axially-stiff cylinder had significantly more 
damage than the other two, i.e., 26.5% at first fiber failure pressure. 

iv. At the first fiber failure pressure it was found that no complete leakage paths developed in any 
of the three cylinders. This fact was particularly remarkable for the case of the axially-stiff cylin- 
der, in spite of the large amount of damage exhibited at that pressure. 
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3.5 Effects of not Degrading C 12 on the Matrix Failure Scenarios 


To answer the question of how the degradation of the inplane shear modulus impacts the 
matrix cracking scenarios for the quasi-isotropic, circumferentially-stiff, and axially-stiff cylin- 
ders, the matrix cracking scenarios obtained by using the baseline degradation scheme are com- 
pared with the resulting matrix failure scenarios when G 12 is not degraded. In both cases the 
maximum stress criterion is used and |3=0.2. 

In Fig. 3.6 the effects of inplane shear modulus degradation on the matrix cracking pattern for 
the quasi-isotropic elliptical cylinder can be seen. The left column of the figure shows the crack- 
ing maps at the pressure levels of 140, 1 70, 200, and 230 psi that correspond to the baseline deg- 
radation scheme shown in Fig. 3.4, which is the scheme reproduced in Fig. 3.5, in the left hand 
column. The right column of Fig. 3.6 shows the cracking maps at the same pressure levels for the 
case of not degrading Gj 2 . As can be observed from the figure, for the quasi-isotropic cylinder, it 
is fair to say that the degradation of the inplane shear modulus has little impact in the matrix 
cracking scenarios at the pressure levels shown, thus concluding that G n is not important for 
resisting matrix crack progression. In fact, at the first three pressure levels there are no noticeable 
differences between the two cases. At 230 psi, though difficult to see, the distribution of the 
cracks is slightly different; however, the overall amount of cracking seems to be the same in both 

cases. For example, in the inner +45° layer, the map for baseline degradation scheme at 230 psi 
shows fewer cracks than the right hand side map where G 12 is not degraded. On the other hand, in 

the inner -45° layer, the former map exhibits more cracks. The argument that the degradation of 
the inplane shear does not play an important role in resisting matrix crack progression for the 
quasi-isotropic cylinder can be further supported by the facts that the first fiber failure pressure 
and the percent damage, listed in Table 3.3, are not affected by the degradation or lack of degrada- 
tion of the inplane shear modulus. 

Figure 3.7 shows the effects of the shear modulus degradation for the circumferentially-stiff 
cylinder. Following a format similar to the one presented for the quasi-isotropic cylinder, the left 
column shows the failure scenarios for the baseline degradation scheme at 140, 170, 200, and 230 
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170 psi 
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Figure 3.6 Effects of shear modulus degradation on the matrix cracking scenario for the 
quasi-isotropic [±45/0/90] s elliptical cylinder, maximum stress criterion, B=0.2 




psi. These cracking maps were previously shown in the middle column of Fig. 3.5. The right col- 
umn shows the cracking maps for the case where G} 2 is not degraded. From the figure, it can be 
seen that for the circumferentially-stiff cylinder, the shear modulus degradation has a small effect 
on the matrix cracking scenarios. Specifically, at all pressure levels shown, the maps for the base- 
line degradation scheme exhibit slightly more cracking than the maps for the case where is 

not degraded. For example, at 140 psi, the left hand side map shows cracks m the 90° layer just 
inside the midsurface, at the sides of the cylinder, s/C &±0. 25. Those cracks do not occur in the 
right hand side map at that pressure. At 1 70 psi, at a glance, there are no significant differences 
between the baseline degradation scheme and the scheme where G ]2 is not degraded. However, 

the innermost 90° layer in the left hand side map shows a few more cracks at s/C&- 0.25. At the 
next pressure level, 200 psi, the same trend occurs, the left hand side map shows more cracks in 

the outer +45° layer near s/C=-0.15 and in the outer and inner -45° layers near s/C=+ 0.25. At 
230 psi pressure level, more cracking still occurs in map for the baseline degradation scheme in 

the ±45° and the inner two 90° layers near the sides of the cylinder. To quantify the impact of 
degrading the inplane shear modulus on the matrix failure scenario for this cylinder; that is, to 
assess to what extent the shear modulus resists the progression of matrix cracks. Table 3.3 shows 
the values for the first fiber failure pressures and the corresponding percent damage associated 
with both cases. As indicated in the table, for the baseline degradation scheme, the percent dam- 
age is slightly larger and the first fiber failure occurs at a slightly lower pressure than when G 12 is 
not degraded. 

With a format similar to the one used for the previous two cylinders. Fig. 3.8 shows the com- 
parison between the matrix failure scenarios for the baseline degradation scheme and the case 
where G )2 is not degraded for the axially-stiff cylinder. The left column corresponding to the 
baseline degradation scheme was previously shown in the furthest right column of Fig. 3.5. 
Though, in Fig. 3.8, at 140, 170, and 200 psi it cannot be seen that the overall amount of cracking 
is larger for the baseline degradation scheme, at 230 psi, the map associated with that scheme 
shows more cracking. At the first two pressure levels, the cracking scenarios are identical for both 
degradation schemes. Interestingly, at 200 psi, the map for the case for the case where G 12 is not 
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degraded shows cracks occurring away from the clamped end of the cylinder in the inner +45° 
layer at s/C»- 0.15. Those cracks do not occur in the left hand side map. Also at this pressure, 

there is more cracking that has progressed axially near the crown in the 0° layer just inside the 
midsurface in the case that neglects the shear modulus degradation. From the maps at 200 psi, it 
would seem that the case where G \ 2 is not degraded exhibits more cracking. However, at 230 psi 
the reverse is true. Longer axial lines of cracks can be identified in the left hand side map in the 

outer +45° layer near the crown, in all the 0° layers near s/C=±0.15, and in the inner +45° layer at 
s/C~-0.15. As with the circumferentially-stiff cylinder, for the axially-stiff cylinder the first fiber 
failure pressure is lower and the associated percent damage is higher for baseline degradation 
scheme than for the case that does not degrade the shear modulus. Those values are indicated in 
Table 3.3. 


Table 3.3: Effects of shear modulus degradation on first fiber failure 
pressures and percent damage, maximum stress criterion, (3=0.2 


Scheme 

baseline degradation 

G \2 not degraded 

Cylinder 

orthotropy 

First fiber 
failure 
pressure, psi 

Percent 

damage 3 

First fiber 
failure 
pressure, psi 

Percent 

damage 3 

quasi-isotropic 

246 

5.5 

246 

5.5 

circumferentially- 

stiff 

253 

6.9 

255 

6.7 

axially stiff 

251 

26.5 

253 

26.2 


a. Number of failed integration points / total number of integration points 
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Figure 3.8 Effects of shear modulus degradation on the matrix cracking scenario for the 
axially-stiff [±45/0 2 /90 1 / 2 ] s elliptical cylinder, maximum stress criterion, p=0.2 





3.6 Influence of (3 on the Matrix Failure Scenarios 

As stated in an earlier chapter, the value of 3=0.2 was adapted from the work in [5]. To inves- 
tigate the effects of the value of the degradation factor on the matrix cracking scenarios for the 
three cylinders studied, the PFA option was configured to use the baseline degradation scheme 
and values of 3=0.05 and 1 .0. It is believed that this range of degradation factors, from severe 
(3=0.05) to no degradation (p=l), would provide a good idea of how sensitive the predicted 
matrix cracking scenarios are to the degradation of the engineering properties. 

Figure 3.9 shows the failure scenarios for the quasi-isotropic cylinder. The first column 
shows the results for the case where the engineering properties are reduced by a factor of 3=0.2. 
These failure maps were last seen in Fig. 3.6. The middle and right columns show the cracking 
maps corresponding to the cases where the engineering properties are reduced by 3=0.05 and 
3=1.0, respectively. As can be seen, if the engineering properties are severely reduced or not 
reduced at all, it makes, practically, no difference in the matrix cracking scenario for the quasi- 
isotropic cylinder. With some difficulty, at 230 psi, it can be observed that fewer cracks are notice- 
able in the furthest right hand side map associated with the case of no degradation. Particularly, in 

layer 5, the 90° layer above the midsurface, near s/C=- 0.25 and in layer 8, a +45° layer, near s/ 
C=+ 0.25. At first fiber failure, the percent damage for the quasi-isotropic cylinder, as indicated in 
Table 3.4, is the same for the three cases. The only effect of using different 3 factors is that the 
first fiber failure pressure seems to increase with less severe degradation. In other words, the fail- 
ure pressure when the engineering properties are not degraded is 2% higher than when the engi- 
neering properties are degraded by 3=0 05. There is no significant difference in the first fiber 
failure pressures between the case of 3=0.2 and the case of 3=0-05. 

In Fig. 3.10, following the same format of Fig. 3.9, the matrix failure scenarios for the circum- 
ferentially-stiff cylinder are shown. Though not easily seen, the failure maps at all pressure levels 
show a slight trend associated with the degree of degradation of the engineering properties. That 
trend, though not pronounced, indicates that the larger the reduction of the engineering properties, 
the more the predicted damage in the cylinder. For example, at 140 psi, the case of 3=0.05 shows 

a few more cracks in the outer +45° layer, layer 9, at s/C~+ 0.25, than the case of 3=0.2. At the 
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Figure 3.9 Effects of degradation factors on the matrix cracking scenario for the quasi-iso- 
tropic [±45/0/90] s elliptical cylinder, maximum stress criterion, baseline degradation scheme 
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same pressure, the map for the case of (3=0.2, in the 90° layer just inside the midsurface, layer 4, 
shows cracks that do not occur in the map for the case of no degradation, or (3=1.0. The same 
trend can be seen in the maps at 1 70 psi. Notice, for the case of (3=0.05, the cracks occurring in the 

inner -45° layer, layer 2, at s/C»+ 0.25. Those cracks do not occur in the map for the case of 

3=0.2. At the same pressure, it can be observed that in the map for 3=0.2, in the outer -45° layer, 
layer 8, at s/C~-0.25 there is more cracking than in the map for the case of no degradation. Simi- 
lar observations can be made for the remainder of the pressure levels. Table 3.4 quantifies that 
trend by indicating the percent damage when first fiber failure occurs. Those values suggest that 
the case of no degradation sustains the least amount of damage, while the case with most severe 
degradation exhibits the most damage. Table 3.4 also indicates that the first fiber failure pressures 
increase with less severe degradation. 

The matrix failure maps for the different degradation factors for the axially-stiff cylinder are 
shown in Fig. 3.11. Again, the format of this figure follows the format of figs. 3.9 and 3. 1 0. Com- 
parison of the maps in Fig. 3.11, shows that as the pressure level increases, the effects of degrad- 
ing the engineering properties by different 3 factors become more evident. Particularly, at 200 and 
230 psi, it is easy to see that the larger the reduction of the engineering properties, the larger the 
damage in the cylinder. At 200 psi, for 3=0.05, more cracks occur, particularly near the crown, 
compared to the map corresponding to 3=0.2. On the other hand, in the case of no degradation, 
there is significantly less cracking than in the case of 3=0.2. A similar scenario can be seen at 230 
psi, where still more cracking near the crown in the case of 3=0-05 can be distinguished com- 
pared to the case of 3=0.2, and significantly less damage is observed in the case of no degrada- 
tion. Table 3.4 lists the first fiber failure pressures and percent damage for the axially-stiff 
cylinder for the three different degradation factors. Those values indicate that, compared to the 
other two cylinders, the degree of degradation of the engineering properties affects the first fiber 
failure pressure and amount of damage for this cylinder. 
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Table 3.4: Effects of degradation factors on first fiber failure pressures and percent damage, 
maximum stress criterion, baseline degradation scheme 


degradation 

factor 

p = 0.05 

P = 0.2 

p = i.o 

Cylinder 

orthotropy 

t->’ m 

riiMiiuci 

failure 

pressure, 

psi 

Percent 

damage 3 

i not iiuwi 

failure 

pressure, 

psi 

Percent 

damage 3 

First fiber 
failure 
pressure, 
psi 

Percent 

damage 3 

quasi-isotropic 

245 

5.5 

246 

5.5 

250 

5.5 

circumferentially- 

stiff 

252 

7.0 

253 

6.9 

258 

6.5 

axially-stiff 

250 

27.3 

251 

26.5 

258 

22.8 


a. Number of failed integration points / total number of integration points 
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3.7 Effects of Recursive Degradation on the Matrix Failure Scenarios 

All the degradation schemes investigated up to this point involved the non-recursive degrada- 
tion of the engineering properties. In a non-recursive degradation, at a given pressure load, the 
engineering properties at the integration points where failure is detected are multiplied by the fac- 
tor P when the failures are encountered. In a recursive degradation, the engineering properties are 
again multiplied by P at each subsequent pressure level. For example, if at the first pressure step, 
failures are detected, the properties at those failed points are multiplied by p. At the second pres- 
sure step, if additional failures are found, the newly failed points are multiplied by P and the pre- 
viously failed points are multiplied by p 2 . To investigate the effects of the recursive degradation 
on the matrix failure scenarios for the cylinders studied, the cracking maps obtained by using the 
PFA option with a recursive degradation of E 2 , G n , v 12 , and v 21 by a factor of p=0.2 are com- 
pared to the results using the baseline degradation scheme with P=0.2. All results are based on the 
maximum stress criterion. It should be noted that use of the recursive scheme results in the degree 
of degradation depending directly on the number of steps taken to reach the final pressure, i.e., 5 
psi increments or 10 psi increments. Wether or not this is an issue could be problem dependent. 

For the quasi-isotropic cylinder. Fig. 3.12 shows the comparisons between the cracking maps 
for the non-recursive and the recursive degradation schemes with P=0.2. Following the format of 
the previous sections, in the left column, the maps corresponding to the non-recursive degradation 
of the engineering properties are shown. The right column correspond to the recursive degrada- 
tion of the properties. It is evident from the figure that the recursive degradation has, for all prac- 
tical purposes, no effect on the failure scenarios for this cylinder. No easily noticeable differences 
can be detected from the maps. At the last two pressure levels, though difficult to see, there are a 
few additional cracks in the recursive case. Particularly, at 200 psi, in the right hand side map, 

there can be seen additional cracking in the innermost 90° layer at s/C»- 0.25. At 230 psi, the 

recursive case shows a few additional cracks in the inner +45° layer near s/C=-0.15. However, 
the differences between the failure maps are insignificant. In fact, the first fiber failure pressures 
and associated percent damage, as indicated in Table 3.5, suggest that this cylinder is not affected 
by the recursive degradation. 
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In Fig. 3.13, the matrix failure scenarios for the non-recursive and recursive degradations for 
the circumferentially-stiff cylinder are shown. From the figure, it is fair to say that, for this cylin- 
der, there are also no significant differences between both schemes; although, in general, the maps 
for the recursive degradation exhibit additional cracking. Specifically, at 170 psi additional crack- 
ing can be seen in the inner -45° layer near the side s/C=+ 0.25 in the map for the recursive case. 
At 200 psi, 2.1s o, 2 few more fnilurcs esn be observed in the map for the recursive ense in purh^u- 

lar, in the inner ±45° layers at both sides near the clamped end. A similar observation can be made 
at 230 psi, where, in layer 3 near the sides at the clamped end, more cracks have progressed cir- 
cumferentially for the recursive degradation. However, quantitatively, the pressure and percent 
damage at first fiber failure listed in Table 3.5 indicates that there is no significant difference 
between the two degradation schemes for the circumferentially-stiff cylinder. 

At a glance, the matrix failure scenarios for the axially-stiff cylinder shown in Fig. 3.14 sug- 
gest that, as for the other two cylinders, there is no significant difference between recursive and 
non-recursive degradation. However, by a close observation, the amount of damage for the axi- 
ally-stiff cylinder due to the recursive degradation is more pronounced than for the other two cyl- 
inders. At 200 psi and 230 psi, several more failures can be observed in the maps for the recursive 

case. At 200 psi for the recursive case, in the crown area in the 0° layers there are more failures 
that have progressed axially. In the right hand side map at 230 psi, there are even more failures 
that have progressed axially and circumferentially. Though difficult to see, this occurs in the top 

two layers near the crown. Note that the gap between the two clusters of failures in the outer ±45° 
layers is smaller for the recursive case, indicating a larger amount of damage in the area surround- 
ing the gap. In the inner +45° layer for the recursive case, it can be seen that there is an additional 
line of cracks in the axial direction near 5/C=-0.05. From the values listed in Table 3.5, it 
becomes clear that for the axially-stiff cylinder, recursive degradation results in more damage at 
first fiber failure and that the pressure at which first fiber failure occurs is lower than for non- 
recursive degradation. 
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Table 3.5: Effects of recursive degradation on first fiber failure pressures 
and percent damage, maximum stress criterion, p=0.2 


Degradation 

Non-recursive 

Recursive 

Cylinder 

orthotropy 

First fiber 
failure 
pressure, psi 

Percent 

damage" 

First fiber 
failure 
pressure, psi 

Percent 

damage* 

quasi-isotropic 

246 

5.5 

245 

5.4 

circumferentially-stiff 

253 

6.9 

252 

7.0 

axially-stiff 

251 

26.5 

249 

27.2 


a. Number of failed integration points / total number of integration points 


For the recursive degradation for all three cylinders, STAGS took considerably longer time to 
run. The number of iterations needed for the analysis to converge was larger than for the non- 
recursive degradation. In particular, for the axially-stiff cylinder, the number of iterations, at the 
higher pressure levels, was above 30, close to the maximum number of iterations allowed by 
STAGS using the default settings for convergence. 
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Figure 3.14 Effect of recursive degradation on the matrix cracking scenario for the axially- 
stiff [±45/0 2 /90|/2l s elliptical cylinder, maximum stress criterion, P=0.2 
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3.8 Evaluation of Computational Algorithms 


A point that can be made at this time is that, in obtaining the results presented here, the PFA in 
STAGS used a specific computational algorithm to compute the stiffness matrix components. 
This algorithm uses the extensional moduli, and E 2 , in a particular fashion in the computation 
of the stiffness matrix components. Specifically, using STAGS notation in the principal material 
coordinate system, the stiffness matrix, which is orthotropic, is given by 


cm{ 1) cm( 4) 0 

cm( 2) cm( 5) 0 

0 0 cm( 9) 


(3.4) 


where the components are computed as follows: 


with 


and 


cm(\) = 


A •£, 


cm( 2) = 


v 12 
A •£, 


cm( 4) = cm(2) 
CM(5) = 
cm( 9) = G 12 


A = 


1 ~ V 12 ' V 21 
E\ • E 2 


(3.5) 


(3.6) 


V 21 = 


(3.7) 


where the superscript ‘o’ indicates that the ratio is based on original, or undegraded, values of the 
extensional moduli. (Note: Equation (3.7) is another way of stating, as was stated in Table 2.1, 
that both Poisson’s ratios are degraded when failure is encountered.) This algorithm can encoun- 
ter numerical singularities if the degradation factor P is very small, or is set to zero to completely 
eliminate either extensional moduli. Although by defining V 2 j using the undegraded ratio helps to 
prevent v 2] from becoming zero or infinity when the degradation factor P is set to zero, the 
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expression for A in Eq. (3.6) becomes very large if either moduli is very small, and infinity if 
either moduli is reduced to zero. In a number of programming languages, including FORTRAN, 
which is the language STAGS utilizes, the very large or infinity value results in a computatio nal 
error, sometimes referred to as overflow. Two options that were recently implemented in STAGS 
circumvent this computational problem. 


The first option, which will be referred to as option 1, defines A without using the extensional 
moduli and computes the matrix components in a different manner than in Eq. (3.5). Specifically, 
this option defines the stiffness matrix components as follows: 


cm{ 1) = 
v 

c/m(4) = — 



cm (2) = cm (4) 


(3.8) 


E 2 

cm{5) = — 
cm{ 9) = G 12 


with 


A = 1 — v, 2 • v 21 (3.9) 

With this approach, setting either extensional moduli to zero presents no particular computational 
problem. 


With the second option, option 2, the components of the undegraded stiffness matrix are mul- 
tiplied by a factor of p. The factor p multiplies the row and column of the stiffness matrix corre- 
sponding to a particular mode of failure. For example, if tensile or compressive matrix failure is 
detected, the elements of the second row and second column of the stiffness matrix are multiplied 
by p. That is, for matrix failure, the stiffness matrix is written as 


cm{\) P • cm(4) 0 

P • cm( 2) P • cm(5) 0 

0 Op- cm{9) 


(3.10) 


If fiber tensile or compressive failure is encountered, the stiffness matrix is then given by 
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(3.11) 


P • cm( 1) p • cm(4 ) 0 

P • cm(2) cm(5 ) 0 

0 OP- cm( 9) 

If shear failure occurs, Eq. (3.4) becomes 


cm(l) cm(4) 0 

cm(2) cm( 5) 0 

0 0 pcm(9) 


(3.12) 


With option 2, the overflow problem is not an issue since the definition of the components cm are 
defined using the undegraded moduli. 


To assess the influence of different computational algorithms on the results corresponding to 
the baseline degradation scheme, with P=0.2, and using the maximum stress criterion, Table 3.6 
shows the number of iterations needed for convergence and percent damage at different pressure 
levels for the three cylinders. As can be seen from the table, no significant differences between 
the three options can be observed for the quasi-isotropic cylinder. In fact, for all practical pur- 
poses, all options go through the same number of iterations and find the same percent damage. 
For the circumferentially-stiff cylinder, the number of iterations varies with algorithm. However, 
it is difficult to determine which option goes through the least number of iterations, since there is 
not a clear trend in one direction or the other. Though the difference is not significant, option 2 
finds the least amount of damage for this cylinder. For the axially-stiff cylinder, the number of 
iterations and percent damage also depend on the algorithm. As for the circumferentially-stiff cyl- 
inder, it is difficult to determine which option goes through the least amount of iterations to con- 
verge, although at the last two pressure levels, option 1 goes through more iterations. It should be 
mention that in running the three options for the axially-stiff cylinder, the extended number of 
iterations (40 maximum) was used. Recall that the default option had to be extended from 22 
maximum iterations to 40 maximum at 190 psi. It should be stated that in the range from 190 to 
240 psi, the extension was needed with all options, since the number of iterations was larger than 
22 at certain pressures. For this cylinder, the difference in percent damage between option 2 and 
the other two options is more evident than the difference found for the circumferentially-stiff cyl- 
inder, particularly, at the last four pressure levels. 
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Figures 3.15, 3.16, and 3.17 show, respectively, the matrix cracking scenarios at different 
pressure levels for the quasi-isotropic, circumferentially-stiff, and axially-stiff cylinders obtained 
using all three computational algorithms. As can be seen from figures 3.15 and 3.16, there are no 
significant differences between the options for the quasi-isotropic and circumferentially-stiff cyl- 
inders. This observation emphasizes the fact, seen in Table 3.6, that all three options predicted 
practically the same amount of damage in these cylinders. In the matrix cracking maps for the axi- 
ally-stiff cylinder, shown in Fig. 3.17, noticeable differences can be observed between option 2 
and the other two options. The maps associated with option 2, at 220 and 240 psi, exhibit the least 
amount of cracking, as also indicated in Table 3.6 by the difference in percent damage. The 
default option and option 1 exhibit no significant differences. Table 3.7 shows the first fiber fail- 
ure pressures and associated percent damage. Interestingly, all options, at the first fiber failure 
pressures, predict basically the same amount of damage for the three cylinders. However, with 
option 2 the first fiber fails at a higher pressure than with the default and option 1. 


Table 3.6: Effects of computational algorithms on the number of iterations and percent 
damage, maximum stress criterion, baseline degradation scheme, (3=0.2 


Cylinder 

orthotropy 

Option 

Default 

1 

2 

Pressure 

Number of 
iterations 

Percent 

damage 3 

Number of 
iterations 

Percent 

damage 3 

Number of 
iterations 

Percent 

damage 3 

quasi-isotopic 

140 

2 

0.0 

2 

0.0 

2 

0.0 

150 

2 

0.0 

3 

0.0 

3 

0.0 

160 

3 

0.1 

3 

0.1 

3 

0.1 

170 

4 

0.2 

3 

0.2 

3 

0.2 

180 

3 

0.4 

3 

0.4 

3 

0.4 

190 

4 

0.6 

4 

0.6 

3 

0.6 

200 

4 

1.0 

4 

1.0 

5 

1.0 

210 

7 

1.8 

7 

1.8 

6 

1.8 

220 

7 

2.8 

6 

2.8 

5 

2.8 

230 

7 

3.9 

6 

3.9 

7 

3.8 

240 

8 

4.9 

7 

4.9 

7 

4.9 
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Table 3.6: Effects of computational algorithms on the number of iterations and percent 
damage, maximum stress criterion, baseline degradation scheme, P=0.2 


Cylinder 

orthotropy 

Option 

Default 

1 

2 

Pressure 

Number of 
iterations 

Percent 

damage 3 

Number of 
iterations 

Percent 

damage 3 

Number of 
iterations 

Percent 

damage 3 

1 

to 

in 

i 

a 

(D 

i— 

o 

o 

i /in 

A ~T V 


n a 

\J 

A 

-T 

n *3 


n *5 

150 

4 

0.6 

5 

0.5 

4 

0.5 


8 

0.9 

8 


7 

0.9 


6 

1.3 

5 

1.3 

7 

1.2 


10 

1.8 

18 

1.8 

6 

1.7 

190 

7 

2.3 

6 

2.4 

6 

2.2 

200 

6 

3.0 

6 

3.0 

6 

2.9 

210 

6 

3.8 

8 

3.8 

8 

3.6 

220 

6 

4.4 

9 

mm 

7 

mm 

230 

8 

5.2 

10 

5.2 

7 

5.0 

240 

8 

5.9 

10 

5.9 

14 

5.8 

M— < 

Cfl 

1 

§ 

140 

2 

0.0 

2 

0.0 

2 

0.0 

150 

2 

0.0 

3 

mom 

2 

0.0 

160 

2 

0.1 

3 


3 

0.1 

170 

3 

0.2 

3 

0.2 

3 

0.2 

180 

3 

0.3 

3 

0.3 

3 

0.3 


24 

0.9 

17 

0.9 

5 

0.9 

200 

18 

4.2 

18 

4.2 

21 

4.2 

210 

22 

10.6 

27 

10.6 

18 

10.0 


220 

29 

15.3 

26 

15.3 

26 

14.7 


230 

19 

19.4 

30 

19.4 

25 

18.4 


240 

25 

22.6 

28 

22.6 

20 

21.5 


a. Number of failed integration points / total number of integration points 
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Table 3.7: Effects of computational algorithms on first fiber failure pressures and 
percent damage, maximum stress criterion, baseline degradation scheme, p=0.2 


Option 

Default 

Option 1 

Option 2 

Cylinder 

First fiber 
failure 
pressure, psi 

Percent 

damage 8 

First fiber 
failure 
pressure, psi 

Percent 

damage 3 

First fiber 
failure 
pressure, psi 

Percent 

damage 3 

quasi- 

isotropic 

246 

5.5 

246 

5.5 

247 

5.5 

circumfere- 

tially-stiff 

253 

6.9 

253 

6.9 

255 

6.8 

axially-stiff 

251 

26.5 

250 

1 26.4 

254 

26.4 


a. Number of failed integration points / total number of integration points 


3.9 Summary of Matrix Failure Scenarios in Internally Pressurized Elliptical 
Cylinders Using Different Degradation Schemes 

In this chapter, the matrix failure scenarios for an elliptical cylinder with quasi-isotropic 
[±45/0/90] 5 , circumferentially-stiff [±45 / 90 2 / 0 l/2 ] s , and axially-stiff [±45 /0 2 / 90 1/2 ] s 

material orthotropies were investigated. The PFA was configured to use the maximum stress fail- 
ure criterion with various degradation schemes. The mode of failure for all cases considered was 
found to be due to excessive tension in the direction perpendicular to the fibers. The pressure and 
locations of the initial matrix failures were found to depend on the material orthotropy. After ini- 
tial failure, as the pressure was increased, it was found that the quasi-isotropic and circumferen- 
tially-stiff cylinders developed similar cracking patterns, although the pressure and the amount of 
damage when first fiber failure occurred was higher for the circumferentially-stiff cylinder. The 
axially-stiff cylinder exhibited a significant amount of damage. The damage developed rapidly 
after about 180 psi. The damage was so dramatic that the default algorithm of the PFA, which 
allows for 22 iterations for convergence, had to be extended to 40 iterations to converge. The per- 
cent damage at first fiber failure was approximately 26.5%. These value was about 5 and 4 times 
larger than for the quasi-isotropic and circumferentially-stiff cylinders, respectively. For all three 
cylinders, it was determined that up to first fiber failure, no leakage paths developed because there 
always remained at least one uncracked layer in the cylinder wall. All these findings agreed well 
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default option 


option 1 


option 2 








Figure 3.15 Matrix cracking scenarios for the quasi-isotropic [±45/0/90] s elliptical cylinder, 
maximum stress criterion, baseline degradation scheme, (3=0.2, different computational 

algorithms 
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option 1 


option 2 






OH 




OH 




OH 



240 psi 


Figure 3.16 Matrix cracking scenarios for the circumferentially-stiff [±45/90 2 /0 1 72 ! s ellipti- 
cal cylinder, maximum stress criterion, baseline degradation scheme, (3=0.2, different com- 
putational algorithms 
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option 1 


option 2 






QH 



QH 




240 psi 


Figure 3.17 Matrix cracking scenarios for the axially-stiff [±45/0 2 /90 1/2 l s elliptical cylinder, 
maximum stress criterion, baseline degradation scheme, P=0.2, different computational 

algorithms 
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with previous works that characterized the behavior and failure of internally pressurized elliptical 
composite cylinders. 

With respect to the importance of the inplane shear modulus degradation, it was found that by 
neglecting the shear modulus in the degradation process there was no difference in the matrix 
cracking scenario for the qua si -isotropic cylinder. For the circumferentially-stiff and axially-stiff 
cylinders, not reducing the shear modulus caused slightly less damage and the first fiber failure 
pressure was found to be slightly higher than the case where the inplane shear modulus was 
degraded. 

! 

The matrix cracking scenarios associated with different degradation factors; specifically, 
P=0.05, P=0.2, and p=l .0, or no degradation, indicated that the larger the reduction of the engi- 
neering properties, the larger the amount of damage. However, the effect was not overly pro- 
nounced, and this trend was not exactly valid for the quasi-isotropic cylinder. The percent damage 
when first fiber failure occurred for the quasi-isotropic cylinder was precisely the same for each 
value of P considered. The pressure at which first fiber failure occurred, however, was highest for 
the case of no degradation and lowest for the lowest value of p. The other two cylinders were 
found to be more sensitive to the value of p. The first fiber failure pressure for the circumferen- 
tially-stiff as well as for the axially-stiff cylinder was found to increase with larger p. The percent 
damage at those pressures actually decreased as p increased. It was interesting to find that for the 
axially-stiff cylinder and the lowest value of P, which led to the largest percent damage at first 
fiber failure, there still remained an uncracked layer, impeding the formation of a leakage path. 

In regard to the effects of recursive degradation of engineering properties, it was found that a 
recursive degradation had an effect similar to the effect of severely reducing the engineering 
properties. This was the case because in a recursive degradation the engineering properties are 
multiplied by the factor P at each load step. For the quasi-isotropic cylinder, it was found that 
there were no significant differences between the cracking maps for a recursive and non-recursive 
degradation. The failure maps for the circumferentially-stiff and axially-stiff cylinders exhibited 
more cracking when a recursive degradation was performed. For the quasi-isotropic and circum- 
ferentially-stiff cylinders, the first fiber failure pressures and percent damage were practically 
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unaffected by the recursive degradation. For the axially-stiff cylinder, the first fiber failure pres- 
sure was found to be lower than the respective pressure for the non-recursive case, and the percent 
damage at first fiber failure was higher for the recursive degradation case. 

Finally, the issue of different computational algorithms within the PFA in STAGS was 
addressed. Because the default algorithm in STAGS may cause a numerical overflow problem if 
the engineering properties are multiplied by a small or zero value of 0, two optional algorithms in 
STAGS may be used to prevent the numerical singularities. The number of iterations for conver- 
gence, the percent damage, and the cracking patterns were found to depend to some degree on the 
computational algorithm used. However, the differences were not significant and one computa- 
tional scheme was not found to be superior. 
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Chapter 4 Matrix Failure Scenarios in Internally 
a rcssarszed Elliptical Cylinders 
Using Different Failure Criteria 


Up to this point, the investigation of the matrix failurd scenarios for the elliptical cylinders 
studied has focused on the use of the maximum stress criterion, which is a non-interactive crite- 
rion. In a non-interactive failure criterion, the conditions of failure only involve one stress compo- 
nent. It is believed that if the matrix failure scenarios are predicted by interactive failure criteria, 
which account for the interaction of different stress components, perhaps more cracking would be 
predicted. If that is the case, a leakage path could be predicted to form before the first fiber fails. 
To investigate if that occurs and assess to what extent the failure predictions depend on the failure 
criterion used, the failure scenarios resulting from using the Hashin and Tsai-Wu failure criteria, 
which are interactive criteria, will be compared to the scenarios provided by the maximum stress 
criterion. 

4.1 Comparison Between the Maximum Stress and Hashin Failure Criteria 

For the Hashin and the maximum stress criteria, the PFA option was set to use the baseline 
degradation scheme with (1=0.2. The maximum stress criterion predictions were discussed previ- 
ously. Interestingly, up to first fiber failure the Hashin criterion predictions were comprised of 
only matrix tensile failures, so all the results presented here involve matrix cracking. The condi- 
tions for matrix tensile failure for the maximum stress and Hashin criteria were defined in Eq. 
(2.6) and Eq. (2.11), respectively. 

Figure 4. 1 shows the matrix failure maps at different pressure levels for the quasi-isotropic 
cylinder. The left column of the figure shows the failure maps corresponding to the maximum 
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stress criterion. These maps were last shown in the left column of Fig. 3.15. The right column dis- 
plays the failure maps predicted by the Hashin criterion. As can be seen, the suspicion that more 
cracking would be predicted by the interactive criterion is confirmed in this case. At 140 psi it can 
be observed that the Hashin criterion predicts two additional failures relative to the predictions of 
the maximum stress criterion. The initial matrix failure pressure according to that criterion occurs 
at approximately 134 psi. This pressure is lower than the pressure predicted by the maximum 
stress criterion, as indicated in Table 4. 1 . The locations of the initial matrix failure, however, are 
predicted to be the same by both criteria. At 170 psi, again more failures are predicted by the 

Hashin criterion, particularly, in the outer ±45° layers. Also, cracks in the inner ±45° layers have 
progressed further circumferentially in the right hand side map. At the next pressure level, 200 
psi, the same trend can be seen, however, the failure distribution predicted by both criteria looks 

very similar. In the right hand side map more cracking occurs in the outer +45° layer at si 

C®+ 0.25, as well as in the 90° layer just inside the midsurface at both sides. The differences 
between the maps, however, are somewhat difficult to detect. At 230 psi, more axial and circum- 
ferential progression of matrix cracks are predicted by the Hashin criterion than by the maximum 

stress criterion. This can be seen in the outer -45° layer near s/C=+0.15, in the inner +45° layer in 

the region of s , /C=-0.15, and in the inner ±45° layers, where more overlapping of cracks occurs 
near the crown. The fact that the Hashin criterion predicts more damage at every pressure level 
than the maximum stress criterion is further confirmed by the percent damage at the first fiber 
failure pressure, as indicated in Table 4.2. Interestingly, both criteria predict the same first fiber 
failure pressure for the quasi-isotropic cylinder. Evidently, though the Hashin criterion predicts a 
higher percentage of damage, the failure distribution is very similar to the one predicted by the 

maximum stress criterion. Consequently, in the Hashin criterion predictions the two 0° layers 
remain uncracked, impeding the formation of a leakage path, as was the situation for the predic- 
tions of the maximum stress criterion. 

Following the same format used in Fig. 4.1, fig 4.2 shows the matrix cracking scenarios pre- 
dicted by both criteria for the circumferentially-stiff cylinder. The failure maps corresponding to 
the maximum stress criterion were last seen in the left column of Fig. 3.16. A close observation of 
Fig. 4.2 reveals that the Hashin criterion also predicts more damage at every pressure level for this 
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cylinder. Interestingly, unlike for the quasi-isotropic cylinder, where the initial matrix failure pres- 
sure predicted by the Hashin criterion was lower than that predicted by the maximum stress crite- 
rion, for this cylinder, both criteria predict the same initial matrix failure pressure, as indicated in 
Table 4.1. Figure 4.2 also shows that for the circumferentially-stiff cylinder both criteria predicts 
slightly different failure distributions. The additional cracks predicted by the Hashin criterion are 
mainly distributed in the circumferential direction. For example, at 140 psi, in the right hand side 

map, more circumferential progression of cracks can be observed in the ±45° layers near si 
C=+0.25. Similar observations can be made for the remaining pressure levels. In particular, at 
230 psi the circumferential progression of cracks predicted by the Hashin criterion is more obvi- 
ous. Note the larger overlapping in the inner ±45° layers and the longer circumferential lines of 
cracks in the inner 90° layer. Regardless the larger amount of cracking and the different failure 

distribution predicted by the Hashin criterion, at 230 psi, the 0° layer can be seen to be uncracked. 
Interestingly, the Hashin criterion predicts that, up to first fiber failure, that layer remains 
uncracked. Therefore, no leakage path forms. 

Figure 4.3 shows the failure predictions for the axially-stiff cylinder. The failure maps pre- 
dicted by the maximum stress criterion were shown previously in the left column of Fig. 3.17. 
Here, again, more cracking can be seen in the failure maps predicted by the Hashin criterion. This 
observation becomes obvious, particularly at 200 and 230 psi, where more cracking can be seen in 
the right hand side maps. At these pressure levels the failure maps display similar failure distribu- 
tions, but a larger amount of cracking occurs in all layers in the right hand side maps. Though it 

cannot be discerned from the map at 230 psi, the 90° layer in the failure scenario predicted by the 
Hashin criterion remains uncracked, regardless of the large damage elsewhere. In fact, up to first 

fiber failure the Hashin criterion predicts that the 90° layer remains crack-free, impeding the for- 
mation of a leakage path. Tables 4.1 and 4.2 list the initial matrix failure pressure and the first 
fiber failure pressure and associated percent damage. 

An interesting observation from the results just discussed is that though the shear stress, x 12 , 
seems to influence the amount of damage predicted by the Hashin criterion, since that stress is 
involved in the matrix tensile failure condition of Eq. (2.11), degrading or not degrading the 
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Figure 4.3 Matrix cracking scenarios for maximum stress and Hashin failure criteria for 
the axially-stiff [±45/0 2 /90 ]/2 ] s elliptical cylinder, baseline degradation scheme, p=0.2 
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Table 4.1: Maximum stress and Hashin failure criteria predictions of 
initial matrix failure pressures 



First matrix failure pressure, psi 

Cylinder 

orthotropy 

Maximum stress 
criterion 

Hashin criterion 

quasi-isotropic 

138 

134 

circumferentially-stiff 

105 

; 105 

axially-stiff 

139 

134 


Table 4.2: Maximum stress and Hashin failure criteria predictions of first fiber failure 
pressures and percent damage, baseline degradation scheme, (5=0.2 


Criterion 

Maximum Stress 

Hashin 

Cylinder 

orthotropy 

First fiber 
failure 
pressure, psi 

Percent 

damage 3 

First fiber 
failure 
pressure, psi 

Percent 

damage 3 

quasi-isotropic 

246 

5.5 

246 

6.2 

circumferentially-stiff 

253 

6.9 

250 

mm 

axially-stiff 

251 

26.5 

252 

30.9 


a. Number of failed integration points / total number of integration points 


inplane shear modulus, G j 2 , as was seen in Figures 3.6, 3.7, and 3.8, did not seem to have a large 
impact when using the maximum stress criterion. 
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4.2 Comparison Between the Maximum Stress and Tsai-Wu Failure Criteria 

As listed in Table 1 , when using the Tsai-Wu failure criterion, the PFA option uses a degrada- 
tion scheme that reduces recursively or non-recursively E 2 , v 12 , and v 21 by a factor of p when 
matrix failure is detected. The shear modulus, G\ 2 , is not reduced. The matrix failure scenarios 
predicted by the maximum stress criterion but not reducing G 12 were previously discussed in sec- 
tion 3.5. Therefore, to assess the differences between the Tsai-Wu criterion and the max imum 
stress criterion, results for that section are compared with the results using the Tsai-Wu criterion. 
A factor of P=0.2 is used. Because the results obtained using the Tsai-Wu criterion predicted only 
matrix tensile failures, again, all the results presented here involve matrix cracking. The condition 
for matrix tensile failure for the Tsai-Wu criterion was defined in Eq. (2.17). 

Figure 4.4 shows the matrix failure scenarios for the quasi-isotropic cylinder predicted using 
both failure criteria. The left column displays the failure maps predicted by the maximum stress 
criterion, while the right column shows the maps predicted by the Tsai-Wu criterion. The left col- 
umn in Fig. 4.4 is the same as the right column of Fig. 3.6. As can be seen from the figure, as it 
was the case for the other interactive criterion, the Tsai-Wu criterion predicts more cracking at all 
pressure levels than the maximum stress criterion. However, unlike the Hashin and the maximum 
stress criteria, the Tsai-Wu criterion predicts the first fiber failure pressure to be lower than 230 
psi. That is why the map at 230 psi in the right column is missing. The first matrix failure pressure 
found by the Tsai-Wu criterion is lower than the pressure found by the maximum stress criterion, 
as indicated in Table 4.3. From Fig. 4.4 it can also be seen that the failure distributions are differ- 
ent for the two criteria; however, the distribution predicted by the Tsai-Wu criterion at 200 psi 
looks similar to the one predicted by the maximum stress at 230 psi. Up to first fiber failure, the 
Tsai-Wu criterion also predicts that two layers remain uncracked. Thus no leakage path forms. 
Table 4.4 shows the first fiber failure pressures and percent damage values predicted by both cri- 
teria. The first fiber failure pressure predicted by the Tsai-Wu criterion is about 12% lower than 
that predicted by the maximum stress criterion. It is interesting to see that when first fiber failure 
occurs, according to the Tsai-Wu criterion, only 3.0% of the cylinder is damaged. This value is 
roughly half the percent damage predicted by the maximum stress criterion. 


85 





The matrix cracking predictions for the circumferentially-stiff cylinder are shown in Fig. 4.5. 
The figure shows that the Tsai-Wu criterion, again, predicts more matrix cracking than the maxi- 
mum stress criterion, which is taken from the left hand column of Fig. 3.16, and a first fiber fail- 
ure pressure lower than 230 psi. As indicated in Table 4.3, the initial matrix failure occurs at a 
similar pressure for both criteria. This was also the case when the Hashin criterion was used. 
Table 4.4 compares the first fiber failure pressures and associated percent damage. As can be seen 
from the table, the first fiber failure pressure predicted by the Tsai-Wu criterion is about 13% less 
than the pressure predicted by the maximum stress criterion. As it was found for the quasi-isotro- 
pic cylinder, the percent damage at first fiber failure predicted by the Tsai-Wu criterion is roughly 
half the value found by the maximum stress criterion. 

Figure 4.6 shows the predictions for the axially-stiff cylinder. The results for the maximum 
stress criterion are taken from the left hand column of Fig. 3. 1 7. For this cylinder, the Tsai-Wu cri- 
terion not only predicts more matrix cracking, but also a different failure distribution. At 140 psi, 
the maps look like the maps, in Fig. 4.4, for the quasi-isotropic cylinder at the same pressure. Also 
the first matrix failure pressures, as indicated in Table 4.3, were similar to the quasi-isotropic 
results. At 170 psi, the failure scenario predicted by the Tsai-Wu criterion exhibits more failures 
and a different failure distribution than the failure scenario predicted by the maximum stress. 

Notice that, in the right hand side map, there can be seen matrix cracks in the inner and outer ±45° 
layers away from the clamped end. At 200 psi, the Tsai-Wu criterion predicts significantly more 
cracks, mainly distributed circumferentially, though there are many cracks in the axial direction. 
The first fiber failure pressures with the associated percent damage are listed in Table 4.4. The 
first fiber failure pressure predicted by the Tsai-Wu criterion is 10% lower than the pressure pre- 
dicted by the maximum stress criterion. There is also a factor of 2 in the percent damage pre- 
dicted. 
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Figure 4.5 Matrix cracking scenarios for maximum stress and Tsai-Wu failure criteria for 
the circumferentially-stiff [±45/90 2 /0 1/2 ] s elliptical cylinder, non-recursive degradation of 

E 2 , vj 2 , and v 21 , P=0.2 


Table 4.3: Maximum stress and Tsai-Wu failure criteria predictions of 
initial matrix failure pressures 



First matrix failure pressure, psi 

Cylinder orthotropy 

Maximum stress 
criterion 

Tsai-Wu 

criterion 

quasi-isotropic 

138 

132 

circumferentially-stiff 

105 

106 

axially-stiff 

139 

133 


Table 4.4: Maximum stress and Tsai-Wu failure criteria predictions of first fiber failure 
pressures and percent damage, non-recursive degradation of f? 2 » v 12 , and V 21» P = 0.2 


Criterion 

Maximum Stress 

Tsai-Wu 

Cylinder 

orthotropy 

First fiber 
failure 
pressure, psi 

Percent 

damage 3 

First fiber 
failure 
pressure, psi 

Percent 

damage 3 

quasi-isotropic 

246 

5.5 

216 

3.0 

circumferentially-stiff 

255 

6.7 

221 

3.9 

axially-stiff 

253 

26.2 

227 

12.2 


a. Number of failed integration points / total number of integration points 


4.3 Summary of Matrix Failure Scenarios in Internally Pressurized Elliptical 
Cylinders Using Different Failure Criteria 

In this chapter, the matrix failure scenarios up to first fiber failure predicted by the non-inter- 
active maximum stress criterion, and two interactive criteria, namely, the Hashin and Tsai-Wu 
failure criteria, were investigated for the quasi-isotropic, circumferentially-stiff, and axially-stiff 
cylinders. The results revealed that all three criteria predicted only failures due to excessive a 2 
stress and similar locations of initial matrix failures. Except for the circumferentially-stiff cylin- 
der, the initial matrix failure pressures, however, were lower when the interactive failure criteria 
were used. For the circumferentially-stiff cylinder, compared to the maximum stress criterion, the 
Hashin failure criterion predicted the same matrix failure pressure and the Tsai-Wu failure crite- 
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non predicted a higher first matrix failure pressure. The failure distributions were found to depend 
on the failure criterion used. The interactive criteria predicted more matrix cracking at all pressure 
levels. Except for the axially-stiff case, the first fiber failure pressures were found to be lower for 
the interactive criteria. Particularly, the first fiber failure pressures were predicted to be signifi- 
cantly lower by the Tsai-Wu criterion. For the axially-stiff case, the Hashin criterion predicted a 
slightly higher first fiber failure pressure. The percent damage at first fiber failure predicted by the 
Hashin criterion was higher than that found by the maximum stress criterion. On the other hand, 
the percent damage at first fiber failure found by the Tsai-Wu criterion was lower than the that 
predicted by the maximum stress criterion. Like the maximum stress criterion, the interactive cri- 
teria predicted that at least one layer remained crack-free, impeding the formation of a leakage 
path. 

As leakage is the primary interest of this study, attention now turns to determining the pres- 
sure levels and failure conditions that lead to leakage. 
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Chapter 5 Failure Analysis Beyond Fiber Failure in 
Internally Pressurized Ellintical Cylinders 

M. %/ 


In the previous chapters, it was shown that regardless the degradation schemes or failure crite- 
ria used to predict the matrix failure scenarios up to first fibler failure, no leakage paths developed 
in any of the cases studied. If the failure analysis continues to higher pressures, at some pressure, 
at certain axial and circumferential locations, there will be failure in some form in all layers, thus 
resulting in a complete leakage path. This chapter will present the results of the failure analysis 
beyond first fiber failure, until a leakage path is found in each of the cylinders. At that condition, 
the pressure, location, and modes of failures associated with leakage are determined and the 
amount of fiber damage assessed. In determining the factors involved with leakage, the PFA in 
STAGS continues to be used for this phase of the study. The PFA configuration employed uses 
only the maximum stress failure criterion with the baseline degradation scheme and a degradation 
factor of P=0.2. It is important to note that in this chapter the definition of the baseline degrada- 
tion scheme also includes the degradation of Fj, Gj 2 , v 12 , and v 21 when fiber failure is encoun- 
tered and G ]2 when inplane shear failure is found. 


5.1 Leakage Path Formation in Internally Pressurized Elliptical Cylinders 

In fig 5.1, the failure progression up to the pressure where a leakage path forms in the quasi- 
isotropic cylinder is shown. The figure is arranged so that the failure maps in the left column are 
associated with matrix failure modes, the middle column corresponds to fiber failure modes, and 
the right column corresponds to the inplane shear mode. The upper left map at a pressure of 230 
psi only displays matrix tensile failures, as shown previously in figs. 3.4, 3.5, and other subse- 
quent figures. The maps at 246 psi correspond to the first fiber failure event. Note that at this pres- 
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sure the map in the middle column shows a single fiber failure, indicated by a blue square symbol, 

at the clamped end, near s/C=-0.15, in the outer +45° layer, layer 8. This failure is due to exces- 
sive compressive aj stress in the fiber direction. This fiber failure is believed to be due to bending 
effects at the clamped end. Also, it should be pointed out that, as mentioned earlier, at the first 

fiber failure pressure, in the matrix failure map, the two 0° layers, layers 3 and 6, are free of fail- 
ures, thus preventing a leakage path formation. At 254 psi, the map in ilie left column indicates 
that the matrix tensile failures continue to progress axially and circumferentially, with new fail- 
ures developing, particularly, in the crown region, away from the clamped end. Those particular 

failures appear in a 0° layer, layer 6, which was previously uncracked. Layer 3 remains 
uncracked. At the same pressure, the map in the middle column shows that more compressive 

fiber failures occur in the outer +45° layer and in layer 6, a 0° layer. Also new tensile fiber fail- 
ures, denoted by the red square symbols, develop in the inner —45° layer, layer 2. The map in the 
right column displays the initial inplane shear failures, denoted by green triangular symbols, 
occurring in layers 3 through 6. It is important to note that these shear failures occur suddenly and 
in the same axial and circumferential region where the initial matrix and fiber failures occurred. 
Although at 254 psi it can be seen that all layers of the laminate sustain failures, there is still not a 
complete leakage path at any particular axial and circumferential location. From the maps at this 
pressure, it becomes evident that a leakage path seems to begin developing at the clamped end 
near s/C=-0.15, where all but one layer, layer 7, exhibit some kind of failure. This can be verified 
if all three maps are superimposed. At the next pressure level shown, 255 psi, the map associated 
with matrix failures displays compressive matrix failures (blue round symbols), due to excessive 

compressive a 2 stress in the direction perpendicular to the fibers, occurring in the +45° layer, 
layer 7, which was the only layer that remained uncracked. The map showing fiber tensile and 
compressive failures indicates that these failures progress circumferentially in the layers where 

they initiated. Also new fiber tensile failures occur in the other 0° layer, layer 3. The inplane shear 
failure map shows that more shear failures develop in most of the layers near j/C=- 0.15. If the 
three maps are superimposed, a complete leakage path can be observed near s/C=-0.15. A combi- 
nation of matrix tensile and compressive, fiber tensile and compressive, and inplane shear failure 
modes are responsible for the formation of this leakage path. The leakage pressure of 255 psi is 
less than 4% greater than the first fiber failure. 
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Figure 5.1 Modes of failure involved in leakage for the quasi-isotropic [±45/0/90] s elliptical 
cylinder, baseline degradation scheme, maximum stress criterion, B=0.2 









For the circumferentially-stiff cylinder, following a similar format to Fig. 5.1, Fig. 5.2 shows 
the failure maps at the pressures of 230 psi, 253 psi, and 254 psi. As shown previously in the cen- 
ter column of Fig. 3.5, at 230 psi, only matrix tensile failures occur. At 253 psi, the initial fiber 
failure occurs. The map in the middle column, at this pressure, shows that, as was the case for the 
\^uuSi~isotrop 10 cylinder, the initial fiber failure occurs due to excessive compressive CTj stress, in 

layer 9, a +45 ° layer, near s/C=-0.15 at the clamped end. Again, bending effects are believed to 
be responsible for this failure in an outer layer. The map in the left column shows a pronounced 
circumferential progression of matrix cracks, particularly, in the three inner layers, relative to the 

230 psi condition. Note also that, in this matrix cracking map, the only 0° layer, layer 5, exhibits 
no damage, preventing the formation of a leakage path. At 254 psi, the initial inplane shear fail- 
ures develop. These shear failures, which appear suddenly, also occur near s/C=-0.15 at the 
clamped end in all layers except in layer 1 . Also at this pressure, more fiber compressive failures 

occur in layer 9 and new fiber tensile failures occur in the inner -45° layer, layer 2. If these three 
maps are superimposed, a complete leakage path can be observed in the region near s/C=-Q. 1 5 at 
the clamped end due to a combination of matrix tensile, fiber compressive and tensile, and inplane 
shear failure modes. The closeness of the first fiber failure pressure, 253 psi, to the leakage pres- 
sure, 254 psi, is quite surprising. 

With the same format of Fig. 5.1 and Fig. 5.2, Fig. 5.3 shows the failure progression that leads 
to the formation of a leakage path in the axially-stiff cylinder. The first pressure level of 230 psi, 
shown previously in the right column of Fig. 3.5, corresponds to a pressure where only matrix ten- 
sile failures occur. At 25 1 psi, as seen in the map of the middle column, initial fiber failure occurs 

at the clamped end near s/C=- 0.15 in the +45° layer, layer 9, due to excessive compressive < 7 j 
stress. As mentioned in previous chapters, though the matrix failure map in the left column exhib- 
its significant amount of damage at first fiber failure pressure, the 90° layer, layer 5, remains 
crack-free, preventing a leakage path formation. The next pressure level shown, 265 psi, is when 
the initial inplane shear failures occur. Again, it can be seen that for this cylinder these failures 
develop suddenly and occur in several layers. More fiber compressive failures and fiber tensile 
failures can also be seen in layers 2, 3, 7, and 9 at 265 psi, in the map of the middle column. Also 
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Figure 5.2 Modes of failure involved in leakage for the circumferentially-stiff 
[±45/90 2 /0 1/2 ] s elliptical cylinder, baseline degradation scheme, 
maximum stress criterion, P=0.2 

at this pressure, though difficult to see in the map of the left column, the single 90° layer remains 
uncracked regardless the large amount of matrix tensile failures. In addition, there is also matrix 

compressive failures occurring in the outer -45° layer, layer 8. Again, if the three maps are super- 
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imposed, a leakage path can be observed at the clamped end near s/C--0.15 due to matrix tensile 
and compressive, fiber tensile and compressive, and inplane shear failure modes. 



• Fiber compressive 
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Inplane shear 
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Figure 5.3 Modes of failure involved in leakage for the axially-stifT [±45/0 2 /90 1/2 ] s elliptical 
cylinder, baseline degradation scheme, maximum stress criterion, (5=0.2 
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5.2 Summary of Failure Analysis Beyond Fiber Failure in Internally 
Pressurized Elliptical Cylinders 


The results presented in this chapter have shown that leakage occurs in the quasi-isotropic, 
circumferentially-stiff, and axially-stiff cylinders studied due to a combination of failure modes. 
Specifically, it was found that matrix tensile and compressive, fiber tensile and compressive, and 
inpiane shear failure modes contribute to the formation of leakage paths. The axial and circumfer- 
ential locations where leakage paths initially form were found to be similar for all three cylinders 
and the leakage pressures were within 4 % of each other. Table 5.1 summarizes the leakage pres- 
sures. This fact suggests that the formation of leakage paljhs depend on the elliptical geometry 
rather than on material orthotropy. 


At leakage, the quasi-isotropic and axially-stiff cylinders sustained similar amount of fiber 
damage. The circumferentially-stiff cylinder, in the other hand, sustained the least amount of fiber 
damage. However, in all the cylinders, when leakage occurs, the overall number of fibers that 
have failed is relatively small. Also, for all cylinders it was found that shear failures develop sud- 
denly, and soon after, leakage occurs. 


Table 5.1: Leakage pressures, baseline degradation scheme, 
maximum stress criterion, (5=0.2 


Cylinder 

orthotropy 

Leakage 

pressure, 

psi 

quasi-isotropic 

255 

circumferentially-stiff 

254 

axially-stiff 

265 
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Chapter 6 Conclusions and Future Work 


6.1 Conclusions 

The goals of the work presented, defined in Chapter 1, have been accomplished. Specifically, 
using a progressive failure analysis approach, the conditions leading to leakage in internally pres- 
surized elliptical cylinders with various material orthotropies have been investigated. In the inves- 
tigation, the progressive failure analysis capability, PFA option, of the finite-element code 
STAGS was utilized. As discussed in Chapter 2, the PFA accounts for the localized loss of stiff- 
ness when material failure occurs by reducing the local engineering properties by a certain degra- 
dation factor. The manner in which this degradation of engineering properties takes place depends 
on the failure modes, which are determined by the use of a failure criterion. The PFA was per- 
formed on a quasi-isotropic, circumferentially-stiff, and axially-stiff cylinders. The results 
obtained using this analysis were organized in two parts. First, the failure scenarios for the three 
cylinders at pressures below the levels to cause fibers to fail were studied. From these failure sce- 
narios obtained using the PFA option with different degradation schemes and various failure cri- 
teria, it was determined that before fibers began to fail in each of the three cylinders, only matrix 
tensile failures, or matrix cracking, develop within the cylinder walls. Though some differences 
between the matrix cracking scenarios were found to depend, in some instances, on the particular 
degradation scheme and failure criterion used, leakage at those pressure levels did not occurred in 
any of the cylinders because at least one layer always remained crack-free, preventing the forma- 
tion of a leakage path. The following list summarizes the main differences determined from the 
comparisons made between the matrix cracking scenarios obtained for the quasi-isotropic, cir- 
cumferentially-stiff, and axially-stiff cylinders using the PFA with several degradation schemes 
and failure criteria. 
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1. PFA with maximum stress failure criterion and several degradation schemes. 

• Effects of material orthotropy: The pressure and locations of initial matrix cracks were 

found to depend on material orthotropy. At higher pressures, the cracking patterns for 
the quasi-isotropic and circumferentially-stiff cylinders were found to be similar. The 
cracking pattern for the axially-stiff cylinder, in the other hand, was remarkably differ- 
ent than the other two cylinders, as the amount of damage sustained by the axially-stiff 
cylinder was significantly higher. 

• Effects of not degrading G ]2 : The importance of the inplane shear modulus in resisting 

matrix cracking was insignificant for the quasi-isotropic cylinder. For the circumferen- 

I 

tially-stifif and axially-stiff cylinders, it was found that by not degrading the shear 
modulus, these cylinders sustain slightly less damage and higher first fiber failure 
pressures. 

• Effects of different degradation factors: A severe degradation of the engineering proper- 

ties was found to produce more damage in the circumferentially-stiff and axially-stiff 
cylinders than a less severe, or no, degradation. This trend was not observed in the 
quasi-isotropic cylinder. In fact, the damage at the first fiber failure pressure for the 
quasi-isotropic cylinder was found to be the same for the different degrees of degrada- 
tion used. 

• Effects of recursive degradation: By continuously degrading the engineering properties 

as the load increased, the cracking scenarios at each load level for the circumferen- 
tially-stiff and axially-stiff cylinders exhibited more damage than when a non-recur- 
sive degradation was performed. However, for the circumferentially-stiff, the amount 
of damage at first fiber failure pressure was, for all purposes, unaffected by the recur- 
sive degradation. For the axially-stiff cylinder, when a recursive degradation was per- 
formed, the first fiber failure occurred at a lower pressure, and the amount of damage 
at this pressure was higher that when a non-recursive degradation was used. The 
matrix failure scenarios for the quasi-isotropic cylinder was unaffected by the recur- 
sive degradation. 

• Effects of different computational algorithms: If a very small or zero value of the degra- 

dation factor is used with the PFA in STAGS, a numerical overflow may occur. Two 
optional algorithms that help prevent that problem were compared with the default 
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algorithm. The failure scenarios for the quasi-isotropic, circumferentially-stiff, and 
axially-stiff cylinders were determined using all three algorithms. No significant dif- 
ferences were found between the different computational algorithms. 

2. PFA with baseline degradation scheme and a degradation factor of 0.2 

• Effects of interactive failure criteria: For the three cylinders studied, the failure scenar- 
ios resulting from using the non-interactive maximum stress criterion were compared 
to the scenarios predicted by the interactive Hashin and Tsai-Wu criteria. Except for 
the circumferentially-stiff cylinder, the initial matrix failure pressures were lower 
when the interactive failure criteria were used. For the circumferentially-stiff cylinder, 
the Tsai-Wu failure criteria predicted a higher first matrix failure pressure. The inter- 
active criteria predicted more matrix cracking at all pressure levels. Except for the cir- 
cumferentially-stiff case, the first fiber failure pressures were found to be lower for the 
interactive criteria. In particular, the first fiber failure pressures were predicted to be 
significantly lower by the Tsai-Wu criterion. For the axially-stiff case, the Hashin cri- 
terion predicted a slightly higher first fiber failure pressure. The amount of damage at 
first fiber failure found by the Tsai-Wu criterion was lower than the that predicted by 
the maximum stress criterion. 

The second part of the results involved the progression of failures occurring at pressures 
above first fiber failure pressures up to leakage pressures. In all three cylinders, it was found that 
leakage occurred due to a combination of matrix, fiber, and inplane shear failure modes in adja- 
cent layers of the laminates. In the circumferentially-stiff and axially-stiff cylinders, at pressures 
when shear failures, which developed suddenly, occurred, complete leakage paths at certain axial 
and circumferential locations formed. In the quasi-isotropic cylinder, soon after shear failures 
appeared, leakage developed. At leakage, the quasi-isotropic and axially-stiff cylinders sustained 
similar amount of fiber damage. The circumferentially-stiff cylinder, in the other hand, sustained 
the least amount of fiber damage. However, in all the cylinders, when leakage occurred, the over- 
all number of fibers that failed was found to be relatively small. The leakage paths formed in 
zones of high bending and high inplane shear stresses, which were believed to cause initial matrix 
failure in the quasi-isotropic and axially stiff cylinders, at the clamped ends of the cylinders. It is 
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suspected the elliptical geometry rather than material orthotropy influenced the formation of these 
paths at those locations. Leakage was found to occur at very similar pressures for the quasi-isotro- 
pic and circumferentially-stiff cylinders. Though the leakage pressure for the axially-stiff cylinder 
was found to be higher, all pressures were within 4% of each other. 

6.2 Future work 

To move forward in this research, the following issues should be addressed in the failure anal- 
ysis of the elliptical cylinder with different material orthotropies: 

I 

i. Since the basic research presented here is based on the idea that non-circular laminated com- 
posite cylinders may be used in future designs of fuel tanks for reusable launch vehicles, a study 
with combined loading, which simulate the conditions under which these tanks may be subjected 
to, could be taken into consideration. 

ii. Also, since typical fuels used in reusable launch vehicles may involve liquefied gases at 
cryogenic temperatures, the basic stages of the research could benefit if thermal stresses are 
included in the analysis. 

iii. Since stress concentrations, which eventually may lead to leakage in non-circular cross 
sections, as shown in this work for the elliptical cylinder, depend on the particular noncircular 
geometry and occur at locations that are not entirely obvious, tailoring of fibers in certain optimal 
fashion could compensate for the non-circular geometry, thus preventing the zones of high 
stresses. 

iv. Experimental work should be undertaken to compare with the results presented here. Criti- 
cal issues to be consider in this regard include not only the fabrication of the elliptical cylinders 
but also the implementation of boundary conditions and the monitoring of failures. 
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